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COUNTING ROOTED TREES: 

THE UNIVERSAL LAW t{n) ~ 

JASON P. BELL, STANLEY N. BURRIS, AND KAREN A. YEATS 


Abstract. Combinatorial classes T that are recursively defined using combinations of the 
standard multiset, sequence, directed cycle and cycle constructions, and their restrictions, 
have generating series T(z) with a positive radius of convergence; for most of these a simple 
test can be used to quickly show that the form of the asymptotics is the same as that for 
the class of rooted trees: , where p is the radius of convergence of T. 


1. Introduction 

The class of rooted trees, perhaps with additional structure as in the planar case, is unique 
among the well studied classes of structures. It is so easy to hud endless possibilities for 
dehning interesting subclasses as the fixpoint of a class construction, where the constructions 
used are combinations of a few standard constructions like sequence, multiset and add-a-root. 
This fortunate situation is based on a simple reconstruction property: removing the root from 
a tree gives a collection of trees (called a forest); and it is trivial to reconstruct the original 
tree from the forest (by adding a root). 

Since we will be frequently referring to rooted trees, and rarely to free (i.e., unrooted) trees, 
from now on we will assume, unless the context says otherwise, that the word ‘tree’ means 
‘rooted tree’. 

1.1. Cayley’s fundamental equation for trees. Cayley [5] initiated the tree investiga¬ 
tions^ in 1857 when he presented the well known inhnite product representation^ 

T(q = 

J>1 

Cayley used this to calculate t{n) for 1 < n < 13. More than a decade later ([7j, [S], pUj 'l 
he used this method to give recursion procedures for finding the coefficients of generating 
functions for the chemical diagrams of certain families of compounds. 
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^This was in the context of an algorithm for expanding partial differential operators. Trees play an 
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^This representation uses t{n) to count the number of trees on n vertices. Cayley actually used t{n) to 
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T{z) = 
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1.2. Polya’s analysis of the generating series for trees. Following on Cayley’s work 
and further contributions by chemists, Polya published his classic 1937 paper^ that presents: 
(1) his group-theoretic approach to enumeration, and (2) the primary analytic technique to 
establish the asymptotics of recursively dehned classes of trees. Let us review the latter as 
it has provided the paradigm for all subsequent investigations into generating series dehned 
by recursion equations. 

Let T(^) be the generating series for the class of all unlabelled trees. Polya hrst converts 
Cayley’s equation to the form 

T{z) = z-exp 

m>l 

From this he quickly deduces that: the radius of convergence p of T( 2 :) is in (0,1) and 
T(p) < cxo. He dehnes the bivariate function 

E(z, w) := ze'^ ■ exp ( ^ T(z™)/m), 

m>2 


giving the recursion equation T = E(z,T). Since E(z,tc) is holomorphic in a neighborhood 
of T he can invoke the Implicit Function Theorem to show that a necessary condition for 
to be a dominant singularity, that is a singularity on the circle of convergence, of T is 

E^(2;,T(^)) = 1. 


From this Polya deduces that T has a unique dominant singularity, namely z = p. Next, 
since ^z{pi T(p)), E^„^(p, T(p)) ^ 0, the Weierstrafi Preparation Theorem shows that p is a 
square-root type singularity. Applying well known results derived from the Cauchy Integral 
Theorem 


( 1 ) 
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■dz 


one has the famous asymptotics 

W tin) 
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which occur so frequently in the study of recursively dehned classes. 


1.3. Subsequent developments. Bender (P, 1974) proposed a general version of the 
Polya result, but Canheld (i, 1983) found a haw in the proof, and proposed a more restricted 
version. Harary, Robinson and Schwenk (im, 1975) gave a 20 step guideline on how to carry 
out a Polya style analysis of a recursion equation. Meir and Moon (|23, 1989) made some 
further proposals on how to modify Bender’s approach; in particular it was found that the 
hypothesis that the coefficients of E be nonnegative was highly desirable, and covered a 
great number of important cases. This nonnegativity condition has continued to hnd favor, 
being used in Odlyzko’s survey paper m and in the forthcoming book US of Flajolet 
and Sedgewick. Odlyzko’s version seems to be a current standard—here it is (with minor 
corrections due to Flajolet and Sedgewick na). 


^Republished in book form in m- 
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Theorem 1 (Odlyzko [22], Theorem 10.6). Suppose 

(2) 'E{z,w) = CijZ^w^ with eoo = 0, eoi < 1, {yi,j) > 0 

i,j>0 

(3) T( 2 ;) = '^^tiz'^ with{'ii)ti>0 

i>l 

are such that 

(a) T(^) is analytic at x = Q 

(b) TW = E(z,T(z)) 

(C) E(2, w) is nonlinear in w 

(d) there are positive integers i,j, k with i < j < k such that 


t j, t j, t/j > 0 

gcd{j - i,k - i) = 1. 


Suppose furthermore that there exist h, r, s > 0 such that 

(e) 'Ei{z,w) is analytic in \z\ < r + 5 and |r(;| < s + 5 

(f) E(r,s) = s 

(g) E^(r,s) = 1 

(h) E;,(r, s) 7 ^ 0 and E^^(r, s) 7 ^ 0 . 

Then r is the radius of convergence ofT, T(r) = s, and as n ^ 00 


rE^{r,s) 

‘luEjyjyjipr ^ s) 




Remark 2. As with Polya’s original result, the asymptotics in these more general theorems 
follow from information gathered on the location and nature of the dominant singularities 
of T. It has become popular to require that the solution T have a unique dominant 
singularity—to guarantee this happens the above theorem has the hypothesis (d). One can 
achieve this with a weaker hypothesis, namely one only needs to require 

(d') gcd {{] -i\i < j and U, tj > 0 }) = 1 . 

Actually, given the other hypotheses of Theorem ^ the condition (d') is necessary and 
sufficient that T have a unique dominant singularity. 


The generalization of Polya’s result that we find most convenient is given in Theorem 1281 
We will also adopt the condition that E have nonnegative coefficients, but point out that 
under this hypothesis the location of the dominant singularities is guite easy to determine. 
Consequently the unique singularity condition is not needed to determine the asymptotics. 

For further remarks on previous variations and generalizations of the work of Polya see §|7| 
The condition that the E have nonnegative coefficients forces us to omit the Set operator 
from our list of standard combinatorial operators. There are a number of complications in 
trying to extend the results of this paper to recursion equations w = G{z,w) where G has 
mixed signs appearing with its coefficients, including the problem of locating the dominant 
singularities of the solution. The situation with mixed signs is discussed in §(21 
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1.4. Goal of this paper. Aside from the proof details that show we do not need to require 
that the solution T have a unique dominant singularity, this paper is not about hnding a 
better way of generalizing Polya’s theorem on trees. Rather the paper is concerned with the 
ubiquity of the form (★) of asymptotics that Polya found for the recursively dehned class of 
trees.^ 

The goal of this paper is to exhibit a very large class of natural and easily recognizable 
operators 0 for which we can guarantee that a solution w = T{z) to the recursion equation 
w = Q{w) has coefficients that satisfy (★). By ‘easily recognizable’ we mean that you only 
have to look at the expression describing 0 —no further analysis is needed. This contrasts 
with the existing literature where one is expected to carry out some calculations to determine 
if the solution T will have certain properties. For example, in Odlyzko’s version. Theorem 
^ there is a great deal of work to be done, starting with checking that the solution T is 
analytic at ^ = 0. 

In the formal specihcation theory for combinatorial classes (see Flajolet and Sedgewick 
ca) one starts with the binary operations of disjoint union and disjoint sum and adds unary 
constructions that transform a collection of objects (like trees) into a collection of objects 
(like forests). Such constructions are admissible if the generating series of the output class 
of the construction is completely determined by the generating series of the input class. 

We want to show that a recursive specihcation using almost any combination of these 
constructions, and others that we will introduce, yield classes whose generating series have 
coefficients that obey the asymptotics (★) of Polya. It is indeed a universal law. The goal 
of this paper is to provide truly practical criteria (Theorem ITHll to verify that many, if not 
most, of the common nonlinear recursion equations lead to (*). Here is a contrived example 
to which this theorem applies: 

(4) w = z + zMSet(Seq{ ^ ^ (2*" + l)(DCyclep,i^, 3 (M;))”. 

nEOdd nEEven 

An easy application of Theorem I7H1 fsee ji l4.29|l tells us this particular recursion equation has 
a recursively dehned solution T{z) with a positive radius of convergence, and the asymptotics 
for the coefficients tn have the form (*). 

The results of this paper apply to any combinatorial situation described by a recursion 
equation of the type studied here. We put our focus on classes of trees because they are by 
far the most popular setting for such equations. 

1.5. First definitions. We start with our basic notation for number systems, power series 
and open discs. 

Definition 3. 

(a) M is the set of reals ; is the set of nonnegative reals. 

(b) P is the set of positive integers. M is the set of nonnegative integers. 

(c) M-‘^[[ 2 ;]] is the set of power series in z with nonnegative coefficients. 

^The motivation for our work came when a colleague, upon seeing the asymptotics of Polya for the first 
time, said “Surely the form (*) hardly ever occurs! (when finding the asymptotics for the solution of an 
equation w = Q{w) that recursively defines a class of trees)”. A quick examination of the literature, a few 
examples, and we were convinced that quite the opposite held, that almost any reasonable class of trees 
defined by a recursive equation that is nonlinear in w would lead to an asymptotic law of Polya’s form (*). 
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(d) Pa is the radius (of convergence) of the power series A . 

(e) For A G M-°[[ 2 ;]] we write A = or A = • 

(f) For r > 0 and zq E C the open disc of radius r about zq is ]D)r(2:o) := {z ■ \z — zo\ < r} 

1.6. Selecting the domain. We want to select a suitable collection of power series to work 
with when determining solutions tc = T of recursion equations w = <h(tc). The intended 
application is that T be a generating series for some collection of combinatorial objects. 
Since generating series have nonnegative coefficients we naturally focus on series in 

There is one restriction that seems most desirable, namely to consider as generating func¬ 
tions only series whose constant term is 0. A generating series T has the coefficient t{n) of 
z"' counting (in some fashion) objects of size n. It has become popular when working with 
combinatorial systems to admit a constant coefficient when it makes a result look simpler, 
for example with permutations we write A(; 2 ) = exp (Q(; 2 )), where A( 2 :) is the exponential 
generating series for permutations, and Q(; 2 ) the exponential generating series for cycles. 
Q( 2 ;) = log(l/(l — z)) will have a constant term 0, but A(; 2 ) = 1/(1 — z) will have the 
constant term 1. Some authors like to introduce an ‘ideal’ object of size 0 to go along with 
this constant term. 

There is a problem with this convention if one wants to look at compositions of operators. 
For example, suppose you wanted to look at sequences of permutations. The natural way to 
write the generating series would be to apply the sequence operator Seq to Ijil — z) above, 
giving X] ~ zY'. Unfortunately this “series” has constant coefficient = cx), so we do 
not have an analytical function. The culprit is the constant 1 in A(; 2 ). If we drop the 1, 
so that we are counting only ‘genuine’ permutations, the generating series for permutations 
is zjiX — z)] applying Seq to this gives zjiX — 2z), an analytical function with radius of 
convergence 1/2. 

Consequently in this paper we return to the older convention of having the constant term 
be 0, so that we are only counting ‘genuine’ objects. 

Definition 4. For A G IR[[2:]] we write A > 0 to say that all coefficients Oi of A are nonneg¬ 
ative. Likewise for B G M[[ 2 :,tc]] we write B > 0 to say all coefficients bij are nonnegative. 
Let 

(a) D0M[2:] := (A G M-°[[ 2 :]] : A(0) = 0}, the set of power series A > 0 with constant 
term 0; and let 

(b) DOM[ 2 :,t(;] := (E G M-°[[ 2 :,tc]] : E(0,0) = 0}, the set of power series E > 0 with 
constant term 0 . Members of this class are called elementary power series.^ 

When working with a member E G DOM[; 2 ,t(;] it will be convenient to use various series 
formats for writing E, namely 

E(;2, w) = CijZ^W^ 
ij 

E(z,w) = 

j 

^We use the name elementary since a recursion equation of the form w = E(z, w) is in the proper form to 
employ the tools of analysis that are presented in the next section. 
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This is permissible from a function-theoretic viewpoint since all coefficients Cij are nonnega¬ 
tive; for any given z,w >0 the three formats converge to the same value (possibly inhnity). 

An immediate advantage of working with series having nonnegative coefficients is that the 
series is defined (possibly inhnite) at its radius of convergence. 

Lemma 5. For T G DOM[z] one has T(pt) ^ [0, cxo]. Suppose T(pt) £ (0, cx)). Then 
Pt < oo ; in particular T is not a polynomial. If furthermore T has integer coefficients then 
Pt < 1- 


2. The theoretical foundations 

We want to show that the series T that are recursively dehned as solutions to functional 
equations w = G{z, w) are such that with remarkably frequency the asymptotics of the coef- 
hcients t^ are given by (★). Our main results deal with the case that G{z, w) is holomorphic 
in a neighborhood of (0,0), and the expansion YhhijZ^w^ is such that all coefficients Pij are 
nonnegative. This covers most of the equations arising from combinations of the popular 
combinatorial operators like Sequence, MultiSet and Cycle. 

The referee noted that we had omitted one popular construction, namely Set, and the 
various restrictions SetM of Set, and asked that we explain this omission. Although the 
equation w = z + 2 Set(tc) has been successfully analyzed in |T7j, there are difficulties when 
one wishes to form composite operators involving Set. These difficulties arise from the fact 
that the resulting equation w = G{z,w) has G with coefficients having mixed signs. A 
general discussion of the mixed signs case is given in 11 1(1.11 and a particular discussion of 
the Set operator in ll l(i.2l Since the issue of mixed signs is so important we introduce the 
following abbreviations. 

Definition 6. A bivariate series E(z,t(;) and the associated functional equation w = E( 2 ;,tc) 
are nonnegative if the coefficients o/E are nonnegative. A bivariate series G{z,w) and the 
associated functional equation w = G{z,w) have mixed signs if some coefficients Pij are 
positive and some are negative. 

To be able to locate the difficulties when working with mixed signs, and to set the stage 
for further research on this topic, we have put together an essentially complete outline of the 
steps we use to prove that a solution T to a functional equation w = E( 2 ;, tc) satisfies the 
Polya asymptotics (*), starting with the bedrock results of analysis such as the Weierstrafi 
Preparation Theorem and the Cauchy Integral Formula. Although this background material 
has often been cited in work on recursive equations, it has never been written down in a single 
unihed comprehensive exposition. Our treatment of this background material goes beyond 
the existing literature to include a precise analysis of the nonnegative recursion equations 
whose solutions have multiple dominant singularities. 

2.1. A method to prove (*). Given E G DOM[z, w] and T G DOM[z] such that T = 
E(z,T), we use the following steps to show that the coefficients tn satisfy (*). 

(a) Show: T has radius of convergence p := px > 0. 

(b) Show: T(p) < oo. 

(c) Show: p < oo. 
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(d) Let: T{z) = where V(0) 7 ^ 0 and gcd{n : v{n) 7 ^ O} = 1. 

(e) Let: u = exp(27ri/g). 

(f) Observe: T{u}z) = u}'^T{z), for |z| < p. 

(g) Show: The set of dominant singnlarities of T is {z ■. z^ = p'^}. 

(h) Show: T satishes a qnadratic eqnation, say 

Qo(2^) + Qi(^)T( 2;) + T(z)^ = 0 

for |; 2 | < p and snfficiently near p, where Qo( 2 ;), Qi( 2 ^) are analytic at p. 

(i) Let: D(z) = Qi(z)2 -4Qo {z), the discriminant oi the eqnation in (g). 

(j) Show: D'(p) 7 ^ 0 in order to conclnde that p is a branch point of order 2, that is, for 
\z\ < p and sufficiently near p one has T( 2 :) = A(p — z) + B(p — z)yjp — z, where A 
and B are analytic at 0, and B(0) < 0. 

(k) Design: A contour that is invariant under multiplication by uj to be used in the 
Cauchy Integral Formula to calculate t{n). 

(l) Show: The full contour integral for t{n) reduces to evaluating the portion lying be¬ 
tween the angles —T^/q and 'njq. 

(m) Optional: One has a Darboux expansion for the asymptotics of t{n). 

Given that E has nonnegative coefficients, items (a)-(f) can be easily established by im¬ 
posing modest conditions on E (see Theorem I28|l . For (g) the method is to show that one has 
a functional equation F(z, T(z)) = 0 holding for \z\ < p and sufficiently near p, that F( 2 ;, w) 
is holomorphic in a neighborhood of (p, T(p)), and that F(p, T(p)) = F^(p, T(p)) = 0, 
but Fi„^(p, T(p)) 7 ^ 0. These hypotheses allow one to apply the WeierstraB Preparation 
Theorem to obtain a quadratic equation for T(z). 

Theorem 7 (WeierstraB Preparation Theorem). Suppose F( 2 ;, ta) is a function of two com¬ 
plex variables and (zo,vJo) is a point in such that: 

(a) F(^,t(;) is holomorphic in a neighborhood of {zq,Wo) 

5F cI^“^F 

(b) F(zo,'«io) = ^(^cwo) = •■■ = —^{zo,wo) = 0 

ow ow'^ ^ 

(c) T—r{Zo,Wo) ^ 0. 

ow'^ 

Then in a neighborhood of {zo,wo) one has F( 2 :, ta) = P( 2 ;, tc)R( 2 :, ta), a product of two 
holomorphic functions P{z, w) andll{z,w) where 

(i) R( 2 :, w) ^ in this neighborhood, 

(ii) P{z,w) is a ‘monic polynomial of degree k’ in w, that is P{z,w) = Qo(^)-l-Qi(- 2 )ta + 
■ ■ ■ + Qk-i{z)w^~^ -|- and the Qi( 2 :) are analytic in a neighborhood of zq. 

Proof. An excellent reference is Markushevich uni, Section 16, p. 105, where one hnds a 
leisurely and complete proof of the WeierstraB Preparation Theorem. □ 

There are two specializations of this result that we will be particularly interested in: k = 1 
gives the Implicit Function Theorem, the best known corollary of the WeierstraB Preparation 
Theorem; and k = 2 gives a quadratic equation for T( 2 ;). 

2.2. k = 1: The implicit function theorem. 

Corollary 8 (k=l: Implicit Function Theorem). Suppose F{z,w) is a function of two com¬ 
plex variables and (zo,uJo) is a point in such that: 
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(a) ¥{z,w) is holomorphic in a neighborhood of {zo,wo) 

(b) F(2 ;o,wo) = 0 

(c) ^(^o,Wo) ^ 0. 

Then there is an e > and a function A( 2 ;) such that for z G 

(i) A(^) is analytic in De(2;o) , 

(ii) F(^,A(2;)) = 0 for z e , 

(iii) for all {z,w) G II 5 £(- 2 o) x De(wo), ifF{z,w) = 0 then w = A(2;). 

Proof. From Theorem [7| there is an e: > 0 and a factorization of 'F{z,w) = L(z, ta)R(z, ta), 
valid in © 2 ( 2 : 0 ) x ©£('W^o), such that R(z, w) 7 ^ 0 for {z, w) G ©^(zo) x ©^(wo), and L(z, w) = 
Lo{z) + w, with Lo( 2:) analytic in ©^(zo)- 

Thus A(;2) = —Lo( 2;) is such that L(^, A(^)) = 0 on ©^(^o); so F(2;,A(;2)) = 0 on 
©(.(zo)- Furthermore, if F(z,t(;) = 0 with {z,w) G ©2(2:0) X ©^(tco), then L(2:,ic) = 0, so 
w = A{z). □ 

2.3. k = 2: The quadratic functional equation. The fact that p is an order 2 branch 
point comes out of the k = 2 case in the WeierstraB Preparation Theorem. 

Corollary 9 (/c = 2). Suppose F( 2 :,tc) is a function of two complex variables and (zo,Wo) 
is a point in such that: 

(a) F(^,t(;) is holomorphic in a neighborhood of {zq,Wq) 

d¥ 

(b) F(2 :o,Wo) = —(2:0,Wo) = 0 
52 F 

(c) ^ 0- 

Then in a neighborhood of (zo,wo) one has F(2:,w) = Q(z, w)R(2:, w), a product of two 
holomorphic functions Q{z,w) and R(2:, w) where 

(i) R( 2 :, w) 7 ^ 0 m this neighborhood, 

(ii) Q(2 :,w) is a ‘monic guadratic polynomial’in w , that is Q{z,w) = Qo( 2:)+ Qi(2:)w + 
w"^, where Qo and Qi are analytic in a neighborhood of Zq. 

2.4. Analyzing the quadratic factor Q,{z,w). Simple calculations are known (see 123 ) 
for hnding all the partial derivatives of Q and R at (zo, 1^0) in terms of the partial derivatives 
of F at the same point. From this we can obtain important information about the coefficients 
of the discriminant D(2:) of Q(z,w). 

Lemma 10. Given the hypotheses (a)-(c) of Corollary let Q{z ,w) and R(z,w) be as 
described in (i)-(ii) of that corollary. Then 

(i) Q(2 :o,Wo) = Q^(2:o,Wo) = 0 

(ii) R(2 :o,Wo) = F^^(2:o,Wo)/2. 

LetY){z) = Qi^zY — 4 :Qq{z) , the discriminant of Q{z,w). 

Then 

(iii) D(2 :o) = 0 

(iv) D'(^o) = -8F2(2 :o,Wo)/F^^(^o,Wo). 
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Proof. For (i) use Corollary (j^I) (b), the fact that R(zo,'W^o) 7^ 0, and 
F(^o,wo) = Q(^o,^i^o)R-(2^o,M^o) 

F^(2;o,m;o) = Q„,(2;o, ^i^o)R-(^o, wo) + Q(2;o, 

= Q^„(zo, w^o)R(2;o,tf^o)- 

For (ii), since Q and vanish and evaluates to 2 at {zo,wo), 

F^^{zo,wo) = 2R(zo,wo)- 

For (hi) we have from (i) 

0 = Qo(^o) + Qi(^o)'M^o + 

0 = Qi(^o) + 2m;o 

and thus 


(5) 

( 6 ) 

From © and we have 
D(zo) 


Qi(2:o) = -2wo 
Qo(2^o) = M^o^- 

Qi(^o)^ ~ 4Qo(^o) = 4:Wo^ — 4wo^ = 0, 


which is claim (iii). 

For claim (iv) start with 

F^(2:o,tf^o) = Q2(2:o,tf^o)R(2:o,tt'o) 

= (Qo(^o) + WoQ'i(2^o))R(^o, Wo). 
From the dehnition of D(z) and ® 

D^(^o) = 2Qi(zo)Qi(^o) — 4Qg(zo) 

= —4(Qg(2;o) + WoQU^o)); 


so 

-4F^(^o,wo) = D'(2 ;o)R(^o, Wo)- 
Now use (ii) to hnish the derivation of (iv). 


□ 


2.5. A square-root continuation of T(^) when z is near p. Let us combine the above 
information into a proposition about a solution to a functional equation. 

Proposition 11. Suppose T G roOM[ 2 ;] is such that 

(a) p:= Pt ^ (0, oo) 

(b) T(p) < cx) 

and F{z,w) is a function of two complex variables such that: 

(c) there is an e > 0 such that F( 2 ;, T( 2 :)) = 0 for \z\ < p and \z — p\ < e 

(d) F(^,t(;) is holomorphic in a neighborhood of (p, T(p)) 

( 9 F 

W F(ftT(p)) =^(p,T(p)) =0 
( 9 F ( 9 ^F 

(f) ^(p.T(p)).^(p.T(p)) > 0. 
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Then there are functions A(z),B(z) analytic at 0 such that 

T( 2 ;) = A(p - z) + B(p - z)^/p - z 
for \z\ < p and near p (see Figure\^, and 


B(0) 


2F4p,T(p)) 

F„„(p,T(p)) 



Figure 1. T{z) = A{p — z) + 'B{p — z)^Jp — z in the shaded region 


Proof. Items (d)-(f) give the the hypotheses of Corollary IHl with ( 2 : 0 , w^o) = (P 5 T(p)). Let 
Qo(^), Qi( 2;) and D(z) = Qi(z)^ — 4Qo(^) be as in Corollary IHl From conclusion (iv) of 
Lemma cni we have 


(7) 


D'(p) = 


g F,(p,T(p)) 

(p, T(p)) 


< 0 . 


From (c) and Corollary IHfi) 

Qo(2^) + Qi(^)T(2;) + T{zY = 0 


holds in a neighborhood oi z = p intersected with Dp(0) (as pictured in Figured]), so in this 
region 

T(p) = -IQi(p) + \PW) 

for a suitable branch of the square root. Expanding D(z) about p gives 


( 8 ) 


D(^) = ^dk{p-zY 


k>l 


since D(p) = 0 by (iii) of Lemma ITUl and di = —D'(p) > 0 by ©. Consequently 

(9) 


T(p) = -iQiW+ ■%/? 


dh 


k>2 


di 


B(p-p) 


holds for \z\ < p and near p. The negative sign of the second term is due to choosing the 
branch of the square root which is consistent with the choice of branch implicit in Lemma 
uni when a = 1/2, given that the t(?7,)’s are nonnegative. 
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Thus we have functions A( 2 ;),B(z) analytic in a neighborhood of 0 with B(0) 7 ^ 0 such 
that 

T( 2 :) = A(p - z) + B(p - z)\/p - z 
for \z\ < p and near p. From diD, m and o 


B(0) = -ivF = -\V-^) 


2F.(p,T(rt) 

F^^(p,T(p)) 


□ 


Now we turn to recursion equations w = E(^, w). So far in our discussion of the role of 
the Weierstrafi Preparation Theorem we have not made any reference to the signs of the 
coefficients in the recursion equation. The following proposition establishes a square-root 
singularity at p, and the proof uses the fact that all coefficients of E are nonnegative. If we 
did not make this assumption then items (USD and (HI below might fail to hold. If (HI 
is false then F^(p, T(p)) may be 0, in which case (*) fails. See section EHl for a further 
discussion of this issue. 


Corollary 12. Suppose T G DOM[z] and E G DOM[. 2 , tc] are such that 

(a) p-.= Pt e (0, 00 ) 

(b) T(p) < cx) 

(c) T(z) = E(z,T(z)) holds as an identity between formal power series, 

(d) E{z ,w) is not linear in w, 

(e) E^ 7^ 0 

(f) {3e > 0) (e(p + £, T(p) + £) < cx)) . 

Then there are functions A(z),B(z) analytic at 0 such that 

T( 2 :) = A{p - z)+ B{p - z)^/p^^ 
for \z\ < p and near p (see Figure\^, and 


B(0) 


2E4p,T(p)) 
Eu,w{p, T(p)) 


Proof. By (f) we can choose e > 0 such that E is holomorphic in 


U — Dp+£( 0 ) X DT(p)+e( 0 ), 


an open 

polydisc neighborhood of 

the graph of T. 

Let 


( 10 ) 

F 

{z,w) := w — 

E{z, 

w). 

Then F 

is holomorphic in U, and one readily sees 

that 


( 11 ) 

F(*,T(z)) = 

T{z) - E{z,T 

(^)) 

= 0 for \z\ < p 

( 12 ) 

F„n,m) = 

1 - E^(^,m;) 



(13) 

F™(p, T(p)) = 

(p, T(p)) 

< 

0 by (d) and E > 0 

(14) 

F.(p,T(p)) = 

-E 4 p,T(p)) 

< 0 

by (e) and E > 0. 


By Pringsheim’s Theorem p is a singularity of T. Thus F^(p, T(p)) = 0 since one cannot 
use the Implicit Function Theorem to analytically continue T at p. 
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We have satisfied the hypotheses of Proposition[TT]— use (IT!^ and (HI to obtain the formula 
for B(0). 


□ 


2.6. Linear recursion equations. In a linear recursion equation 

w = Ao(z) + Ai( 2 :)ta 


one has 
(15) 


w 


Ao(^) 

1 - Ai(z)' 


From this we see that the collection of solutions to linear equations covers an enormous 
range. For example, in the case 

w = Ao( 2 ;) + zw, 

any T{z) G DOM[z] with nondecreasing eventually positive coefficients is a solution to the 
above linear equation (which satisfies Ao( 2 ;) + zw > Q) if we choose Ao(^) := (1 — z)T{z). 

When one moves to a Q{w) that is nonlinear in w, the range of solutions seems to be 
greatly constricted. In particular with remarkable frequency one encounters solutions T(z) 
whose coefficients are asymptotic to C. 


2.7. Binomial coefficients. The asymptotics for the coefficients in the binomial expansion 
of {p — z)" are the ultimate basis for the universal law (*■). Of course if a G N then {p — z)" 
is just a polynomial and the coefficients are eventually 0 . 


Lemma 13 (See Wilf [211, P- 179). For a G M \ N and p G (0, cx)) 

M(p-z)" = (-1) 

2.8. The Flajolet and Odlyzko singularity analysis. In HU Flajolet and Odlyzko de¬ 
velop frans/er theorems via singularity analysis for functions 8 ( 2 ;) that have a unique dom¬ 
inant singularity. The goal is to develop a catalog of translations, or transfers, that say: if 
8 ( 2 ;) behaves like such and such near the singularity p then the coefficients s{n) have such 
and such asymptotic behaviour. 

Their work is based on applying the Cauchy Integral Formula to an analytic continuation of 
8 (^) beyond its circle of convergence. This leads to their basic notion of a Delta neighborhood 
A of p, that is, a closed disc which is somewhat larger than the disc of radius p, but with 
an open pie shaped wedge cut out at the point z = p (see Fig. El). We are particularly 
interested in their transfer theorem that directly generalizes the binomial asymptotics given 
in Lemma 


a 


n 


^a—n 
P ~ 




F(-«) 


p-^n-^-\ 


Proposition 14 (Hll) Corollary 2). Let p G (0, cx)) and suppose 8 is analytic in A \ {p} 
where A is a Delta neighborhood of p. If a ^ N and 

(16) 8 ( 2 ;) ~ K(^p — z)^ 


as z ^ p in A, then 

s{n) [z'^] K(^p — z)°‘ 







Kp^ 

r(-«) 


■p-^n 


—a—1 
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Delta region for a single 
singularity 



The corresponding 
contour shape 


Figure 2. A Delta region and associated contour 


Let us apply this to the square-root singularities that we are working with to see that one 
ends up with the asymptotics satisfying (*). 


Corollary 15. Suppose S G DOM[ 2 ;] has radius of convergence p G (0, cx)), and p is the only 
dominant singularity of S. Furthermore suppose A and B are analytic at 0 with B(0) < 0, 
A(0) > 0 and 

(17) S{z) = A{p - z) + B(p - z)\/p- z 

for z in some neighborhood of p, and \z\ < p. 

Then 

-B(0)v/p „ 3/2 

20F ■ 


s{n) ~ [z”] B(0)x/p- z 


Proof. One can find a Delta neighborhood A of p (as in Fig. ED such that S has an analytic 
continuation to A \ {p}; and for z G A and near p one has (ITTp holding. Consequently 

S(z)-A(0) ~ 

as p in A. This means we can apply Proposition El to obtain 


s(n) 


- 3/2 

r(-i/2) P 


□ 


2.9. On the condition B(0) < 0. In the previous corollary suppose that B(0) = 0 but 
B 7 ^ 0. Let bk be the first nonzero coefficient of B. The asymptotics for s(?7,) are 


s{n) ~ hk[z^]{p-z) 




giving a law of the form Cp~"‘n~^~ 2 . We do not know of an example of S defined by a 
nonlinear functional eguation that gives rise to such a solution with k > Q, that is, with 
the exponent of n being —^/2, or —7/2, etc. Meir and Moon (p. 83 of |2I], 1989) give the 
example 


w = 


{l/Q)e'^^z^/n'^ 


n>l 


where the solution w = T has coefficient asymptotics given by tn ^ C/n. 
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2.10. Handling multiple dominant singularities. We want to generalize Proposition 
d to cover the case of several dominant singnlarities eqnally spaced aronnd the circle of 
convergence and with the fnnction S enjoying a certain kind of symmetry. 

Proposition 16. Given q and p G (0, cx)) let 

uj := 

Uq,p := {(nV: j = 0, 

Suppose A is a generalized Delta-neighborhood of p with wedges removed at the points in Uq^p 
(see Fig. O/or q = 3), suppose S is continuous on A and analytic in A \ Uq^p, and suppose 




A Delta region for The corresponding 

3 singularities contour shape 


Figure 3. Multiple dominant singnlarities 


d is a nonnegative integer such that S[u!z) = a;‘^S( 2 ;) for z E A. 
IfS{z) rsj K{p — z)" as z ^ p in A and a ^ N then 

qKp'^ 


s{n) 


r(-a) 




if n = d mod q, 


s{n) = 0 otherwise. 

Proof. Given £ > 0 choose the contour C to follow the boundary of A except for a radius £ 
circular detour around each singularity uj^p (see Fig. |2I). Then 

s(n) = J- f ^dz. 


27ri 


n+1 


Subdivide C into q congruent pieces Cq, ..., Cg_i with Cj centered around cn-’p, choosing as 
the dividing points on C the bisecting points between successive singularities (see Fig. lUfor 
q = 3). Then Cj = u^Cq. Let Sj(n) be the portion of the integral for s(n) taken over Cj, that 
is: 


Sj[n) = 


27ri Jc, z'^+^ 


dz. 


Then from S{ujz) = a;‘^S(z) and Cj = uj^Cq we have 

- f S(z) 


Sj[n) = 


27ri z^~^^ 


dz 


1 

27ri 


'Co 
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uj'^^S{z) 


ui^dz 













Figure 4. The congruent contour segments Cj 


so 

s(n) = 


We have reduced the integral calculation to the integral over Co, and this proceeds exactly 
as in in the unique singularity case described in Propositional D 

Let us apply this result to the case of S(z) having multiple dominant singularities, equally 
spaced on the circle of convergence, with a square-root singularity at p. 

Corollary 17. Given g G P and p G (0, cx)) let 

w := 

Uq,p := j = 0, 

Suppose S G D0M[2;] has radius of convergence p G (0, cxo), Uq^p is the set of dominant 
singularities of S, and S(a; 2 ;) = a;‘^S(^) for \z\ < p and for some d eN. 

Furthermore suppose A and B are analytic at 0 with B(0) < 0, A(0) > 0 and 

(18) S(z) = A(p - z) + B(p - z)^Jp - z 
for z in some neighborhood of p, and \z\ < p. Then 

(19) s{n) ~ ^ n = d mod q. 

Otherwise s{n) = 0. 


= u:F<i-n)±_ / 

2vrz Jcq 


q-l 

^sj(n) 
j=0 
q-l 

^u;^(^--))so(n) 
j=0 

qso(n) if n = d mod q 
0 otherwise. 
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Proof. Since the set of dominant singularities Uq^p is finite one can find a generalized Delta 
neighborhood A of p (as in Fig. ED such that S has a continuous extension to A which 
is an analytic continuation to A \ Uq^p] and for z E A and near p one has m holding. 
Consequently 

S(z)-A(0) ~ B(0)^/p^ 

as z —p in A. This means we can apply Propositionto obtain dTHl) . □ 

2.11. Darboux’s expansion. In 1878 Darboux published a procedure for expressing 
the asymptotics of the coefficients s{n) of a power series S with algebraic dominant singular¬ 
ities. Let us focus hrst on the case that S has a single dominant singularity, namely z = p, 
and it is of square-root type, say 

S( 2 r) = A{p - z)+B{p - z)^/p^^ 

for \z\ < p and sufficiently close to p, where A and B are analytic at 0 and B(0) < 0. From 
Proposition ITTI we know that 

s(n) = {l + o{l))b{f))[z'^]^/p^. 

Rewriting the expression for S(z) as 

OO 

S(^) = '^{aj{p-zY + hj{p-zy^'^^ 
i=o 

we can see that the mth derivative of S ‘blows up’ as z approaches p because the mth 
derivative of the terms on the right with j < m involve terms with p — z to a negative power. 
However for j > m the terms on the right have mth derivatives that behave nicely near p. 
By shifting the troublesome terms to the left side of the equation, giving 

S^(z) := S(z) - (aj{p - zy + bj{p - zy+^^ 

j<m 

= Y > 

j>m 

one can see by looking at the right side that the mth derivative of Sm(^) has a 

square-root singularity at p provided some bj 7 ^ 0 for j > m. Indeed (z) is very much 
like 8 ( 2 ;), being analytic for \z\ < p provided z y p. If 6 ^ 7 ^ 0 we can apply Proposition EU 
to obtain (for suitable Cm) 

7”]SSrP) ~ C„p-”n-i 

and thus 

[z^]Smiz) ~ Cmp-^n-^-l 

This tells us that 

(«j(p-2:)^ + &i(p-+ (1+ o(l))Cmp“”n“'""t. 

j<m 

For n > m the part with the Oj drops out, so we have the Darboux expansion 





The case of multiple dominant singularities is handled as previously. Here is the result for 
the general exponent a. 


Proposition 18 (Multi Singularity Darboux Expansion). Given g G P /ef 

uj := 

Uq,p := j = 0, 

Suppose we have a generalized Delta-neighborhood A with wedges removed at the points in 
Uq^p (see Fig. 01 and S is analytic in A \ Uq^p. Furthermore suppose d is a nonnegative 
integer such that S(a;z) = a;^S(z) for \z\ < p. 

If 

S{z) = A{p-z) + B{p-z){p-zr 

for \z\ < p and in a neighborhood of p, and a ^ N, then given m G N with bm 0 there is a 
Cm 7 ^ 0 such that for n = d mod q 

s{n) = + 

j<m 


2.12. An alternative approach: reduction to the aperiodic case. In the literature 
one hnds references to the option of using the aperiodic reduction V of T, that is, using 
T{z) = z^Y{z^) where V(0) 7 ^ 0 and gcd{n : v{n) 7 ^ 0} = 1. V has a unique dominant 
singularity at pv = Pt^-, so the hope would be that one could use a well known result like 
Theorem 0 to prove that (★) holds for v{n). Then t{nq + d) = v{n) gives the asymptotics 
for the coefficients of T. 

One can indeed make the transition from T = E(z, T) to a functional equation V = 
H{z,V), but it is not clear if the property that E is holomorphic at the endpoint of the 
graph of T implies H is holomorphic at the endpoint of the graph of V. Instead of the 
property 

(3£ > 0) ^E(p + E, T(p) + £:) < 00 ^ 
of E used previously, a stronger version seems to be needed, namely: 


(Vp > 0) E(p,p) < cxo (3e > 0) (e(p + £,p + e) < cxo 


We chose the singularity analysis because it sufficed to require the weaker condition that 
E be holomorphic at (p, T(p)), and because the expression for the constant term in the 
asymptotics was far simpler that what we obtained through the use of V = H(z, V). Fur¬ 
thermore, in any attempt to extend the analysis of the asymptotics to other cases of recursion 
of equations one would like to have the ultimate foundations of the WeierstraB Preparation 
Theorem and the Cauchy Integral Theorem to fall back on. 


3. The Dominant Singularities of T{z) 

The recursion equations w = B{z, w) we consider will be such that the solution w = T has 
a radius of convergence p in ( 0 , cx)) and finitely many dominant singularities, that is hnitely 
many singularities on the circle of convergence. In such cases the primary technique to hnd 
the asymptotics for the coefficients t(n) is to apply Cauchy’s Integral Theorem (0). Experi¬ 
ence suggests that properly designed contours C will concentrate the value of the integral (0) 
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on small portions of the contour near the dominant singularities of T—consequently great 
value is placed on locating the dominant singularities of T. 

Definition 19. For T G ID)OM[ 2 ;] with radius p G (0, cx)) let DomSing(T) he the set of 
dominant singularities of T, that is, the set of singularities on the circle of convergence of 

T. 

3.1. The spectrnm of a power series. 

Definition 20. For A G DOM[^] let the spectrum Spec(A) of A be the set of n such that 
the nth coefficient a{n) is not zero.^ It will be convenient to denote Spec(A) simply by A, so 
we have 

A = Spec(A) = {n : a{n) ^ 0}. 

In our analysis of the dominant singularities of T it will be most convenient to have a 
simple calculus to work with the spectra of power series. 

3.2. An algebra of sets. The spectrum of a power series from ]D)OM[z] is a subset of positive 
integers; the calculus we use has certain operations on the subsets of the nonnegative integers. 

Definition 21. For /, J C N and j, m G N let 

I + J 
1-3 
m ■ J 
0 © J 
m Q J 

IQJ 
m\J 

3.3. The periodicity constants. Periodicity plays an important role in determining the 
dominant singularities. For example the generating series '~F{z) of planar (0,2)-binary trees, 
that is, planar trees where each node has 0 or 2 successors, is dehned by 

T{z) = z + zT{z)\ 

It is clear that all such trees have odd size, so one has 

CO OO 

T(^) = = z'^t{2j + l){zy. 

j=0 j=0 

®In the 1950s the logician Scholz defined the spectrum of a first-order sentence y to be the set of sizes of 
the finite models of y. For example if y is an axiom for fields, then the spectrum would be the set of powers 
of primes. There are many papers on this topic: a famous open problem due to Asser asks if the collection of 
spectra of first-order sentences is closed under complementation. This turns out to be equivalent to an open 
question in complexity theory. The recent paper m of Fischer and Makowsky has an excellent bibliography 
of 62 items on the subject of spectra. 

For our purposes, if A( 2 ;) is a generating series for a class A of combinatorial objects then the set of sizes 
of the objects in A is precisely Spec(A). 


:= {i + j :iel,j eJ} 

:= — j : i G /, } where j < min(J) 

:= {m ■ j : j E J} for m > 1 
■■= { 0 } 

:= J + ■ y + for m > 1 

m—times 

:= IJfOJ 

i&I 

^ (Vj G J) (m|j). 


18 





This says we can write T(z) in the form 

T{z) = zV(z^). 

From snch considerations one hnds that T(z) has exactly two dominant singnlarities, p and 
—p. (The general result is given in Lemma 1221) 

Lemma 22. For A G DOM[ 2 ;] let 

p := gcdA d := min A q := gcd(A —d). 

Then there are \d{z) and V( 2 ;) in R-°[[^]] such that 

(a) A(^) = \j[zP) with gcd{U) = 1 

(b) aIz) = z'^V{z'^) with V(0) ^ 0 and gcd(l/) = 1. 

Proof. (Straightforward.) □ 

Definition 23. With the notation of Lemma\EE U(^^) is the purely periodic form of A{z); 
and z'^V ( 2 ;'?) is the shift periodic form of A(z). 

The next lemma is quite important—it says that the q equally spaced points on the circle 
of convergence are all dominant singularities of T. Our main results depend heavily on the 
fact that the equations we consider are such that these are the only dominant singularities 
of T. 

Lemma 24. Let T G DOM[z] have radius of convergence p G (0, cxo) and the shift periodic 
formz'^\{z'^). Then 

{z : z'^ = p'^} C DomSing(T). 

Proof. Suppose zq‘^ = p‘^ and suppose S(z) is an analytic continuation of T( 2 ;) into a neigh¬ 
borhood De(2;o) of zq. Let oj := zo/p. Then o;^ = 1. The function So(;2) := S(a;^)/a;‘^ is an 
analytic function on D£(p). For ^ G De(p) fl Dp(0) we have 

u!Z G D£(zo) n Dp(0), 

so 

So(^) = S(a;z)/a;'^ = T{ujz)/u‘^ = T{z). 

This means So(;s) is an analytic continuation of T{z) z = p, contradicting Pringsheim’s 
Theorem that p is a dominant singularity. □ 

3.4. Determining the shift periodic parameters from E. 

Lemma 25. Suppose T(z) = E( 2 :, T( 2 :)) is a formal recursion that defines T G DOM[z], 
where E GDOM[z,t(;]. Let the shift periodic form of T{z) be z‘^W{z'^). Then 

d = min(T) = min(Eo) 

q = gcd(T-d) = gcd IJ -f (n - l)dj. 

n>0 

Proof. Since T is recursively dehned by 

T(^) = E,.(^)T(2)” 

n>0 
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one has the hrst nonzero coefficient of T being the hrst nonzero coefficient of Eq, and thus 
d = niin(T) = niin(Eo). It is easy to see that we also have q = gcd(T — d). 

Next apply the spectrum operator to the above functional equation to obtain the set 
equation 

T = y + n © T, 

n>0 

and thus 

T-d = (^En +{n-l)d + nQ {T - d)^. 

n>0 

Since q = gcd(T — d) it follows that q\r := gcd ^ IJ^ + (n — l)(ij. 

To show that r\q, and hence that r = q, note that 

tc = (^En + {n - l)d + n Q v?j 

n>0 

is a recursion equation whose unique solution is w = T — d. Furthermore we can hnd the 
solution w by iteratively applying the set operator 

0(tc) := (^En +{n-l)d + nQv?j 

n>0 

to 0, that is, 

T-d = lim 0'^(0). 

n—»'00 

Clearly r | 0, and a simple induction shows that for every n we have r | 0”(0). Thus r | (T — 
d), so r I g, giving r = q. This hnishes the proof that q is the gcd of the set IJ^ + (n — 

l)d). □ 

3.5. Determination of the dominant singularities. The following lemma completely 
determines the dominant singularities of T. 

Lemma 26. Suppose 

(a) T G DOM[^] has radius of convergence p G (0, cxo) with T(p) < oo, and 

(b) T(z) = E(z,T(; 2 )), where E G roOM[z,tc] is nonlinear in w and holomorphic on (the 
graph of)T. 

Let the shift periodic form of'T{z) he Then 

DomSing(T) = {z ■. z^^ = 

Proof. By the usual application of the implicit function theorem, if 2 ; is a dominant singularity 
of T then 

(20) E^(z,T(z)) =1. 

As p is a dominant singularity we can replace (OHll by 

(21) E^(z,T( 2;)) = E^(p,T(p)). 
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Let U{zP) be the purely periodic form of T(z)). As the coefficients of E^,, 

nonnegative it follows that m implies 

DomSing(T) C [z : = //*}. 

We know from Lemma 1^ that 

{z : z'^ = C DomSing(T), 

consequently q\p. 

To show that p < q hrst note that if m G N then 

gcd(m + T)\q. 

For if r = gcd(m + T) then for any n G T we have r|(m + n) and r|(m + d). Consequently 
r|(n — d), so r|(T — d), and thus r\q. 

Since 

U(z^) = E^{z,T{z)) = J2^n{z)nT{zr-\ 

n>l 

applying the spectrum operator gives 

Spec(U(z^’)) = ljE„ + (n-l)0T. 

n>l 


Choose n >2 such that En ^ 0 and choose a E En. Then 

Spec(U(z^’)) 0 E„ + (n-l)©T 
0 (a + (n — 2)d) + T, 

so taking the gcd of both sides gives 

p = gcd Spec (U( 2 ;^)) 


< 


gcd ^(a + (n - 2)d) + 




With p = q it follows that we have proved the dominant singularities are as claimed. 


□ 


3.6. Solutions that converge at the radius of convergence. The equations w = 0{w) 
that we are pursuing will have a solution T that converges at the hnite and positive radius 
of convergence px- 


Definition 27. Let 

DOM*[^] := {T G D0M[2;] : px G (0, cxd), T(px) < cxd}. 

3.7. A basic theorem. The next theorem summarizes what we need from the preceding 
discussions to show that T = E(z, T) leads to (★) holding for the coefficients of T. 

Theorem 28. Suppose T G ©OM[^] and E G DOM[^, tc] are such that 

(a) T(a) = EpTW) holds as an identity between formal power series 

(b) T G ©OM*[z] 

(c) E{z, w) is nonlinear in w 

(d) E, ^ 0 

(e) (3£ > 0) ^E(p + e, T(p) + e) < cx)j . 
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Then 


t{n) 


pE,{p,T{p)) 


2vrE^^(p, T(p)) 
Otherwise t{n) = 0. Thus (*) holds on {n : t{n) > 0} 
Proof. By Corollary [T21 Corollary HU and Lemma M 


p ""n for n = d mod q. 


□ 


4. Recursion Equations using Operators 

Throughout the theoretical section, §12 we only considered recursive equations based on 
elementary operators E{z,w). Now we want to expand beyond these to include recursions 
that are based on popular combinatorial constructions used with classes of unlabelled struc¬ 
tures. As an umbrella concept to create these various recursions we introduce the notion of 
operators 0. 

Actually if one is only interested in working with classes of labelled structures then it 
seems that the recursive equations based on elementary power series are all that one needs. 
However, when working with classes of unlabelled structures, the natural way of writing down 
an equation corresponding to a recursive specihcation is in terms of combinatorial operators 
like MSet and Seq. The resulting equation w = 0(tc), if properly designed, will have a 
unique solution T(z) whose coefficients are recursively dehned, and this solution will likely 
be needed to construct the translation oiw = Q{w) to an elementary recursion w = E(^, w), 
a translation that is needed in order to apply the theoretical machinery of §12 

4.1. Operators. The mappings on generating series corresponding to combinatorial con¬ 
structions are called operators. But we want to go beyond the obvious and include complex 
combinations of elementary and combinatorial operators. For this purpose we introduce a 
very general dehnition of an operator. 

Definition 29. An operator is a mapping 0 : DOM[ 2 ;] —>• DOM[ 2 ;] . 

Note that operators 0 act on DOM[z] , the set of formal power series with nonnegative 
coefficients and constant term 0. As mentioned before, the constraint that the constant terms 
of the power series be 0 makes for an elegant theory because compositions of operators are 
always dehned. 

A primary concern, as in the original work of Polya, is to be able to handle combinatorial 
operators 0 that, when acting on T(z), introduce terms like T{z‘^),T{z^) etc. For such 
operators it is natural to use power series T(z) with integer coefficients as one is usually 
working in the context of ordinary generating functions. In such cases one has p < 1 for the 
radius of convergence of T, provided T is not a polynomial. 

Definition 30. An integral operator is a mapping 0 : IDOM[ 2 ;] —> roOM[ 2 :] , where 
IDOM[z] := |A G N[[2:]] : A(0) = O}, the set of power series with nonnegative inte¬ 
ger coefficients and constant term zero. 

Remark 31. Many of the lemmas, etc, that follow have both a version for general operators 
and a version for integral operators. We will usually just state and prove the general version, 
leaving the completely parallel integral version as a routine exercise. 
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4.2. The arithmetical operations on operators. The operations of addition, multipli¬ 
cation, positive scalar multiplication and composition are defined on the set of operators in 
the natural manner: 

Definition 32. 

(01 + 02)(T) := 0i(T) + 02(T) 

(01 ■ 02)(T) := 0i(T) ■ 02(T) 

(c-0)(T) := c-0(T) 

(0io02)(T) := 0i(02(T)), 

where the operations on the right side are the operations of formal power series. A set of 
operators is closed if it is closed under the four arithmetical operations. 

Note that when working with integral operators the scalars should be positive integers. 
The operation of addition corresponds to the construction disjoint union and the operation 
of product to the construction disjoint sum, for both the unlabelled and the labelled case. 
Clearly the set of all [integral] operators is closed. 

4.3. Elementary operators. In a most natural way we can think of elementary power 
series E( 2 :,t(;) as operators. 

Definition 33. Given E( 2 :,t(;) G DOM[ 2 :, tc] let the associated elementary operator be given 
by 

E : T E{z, T) for T e DOM. 

Two particular kinds of elementary operators are as follows. 

Definition 34. Let A G DOM[;2]. 

(a) The constant operator 0a is given by 0a ; T i—>• A for T G DOM[; 2 ], and 

(b) the simple operator A{w) maps T G DOM[; 2 ] to the power series that is the formal 
expansion of 

n>l j>l 

4.4. Open elementary operators. 

Definition 35. Given a,6 > 0, an elementary operator E(; 2 ,t(;) is open at {a,b) if 

(3e > 0) ^E(a + e,b + e) < cx)j. 

E is open if it is open at any a,b > 0 for which E(a, b) < oo. 

Eventually we will be wanting an elementary operator to be open at (p, T(p)]) in order to 
invoke the WeierstraB Preparation Theorem. 

Lemma 36. Suppose A G DOM[ 2 :] and a,b > 0. 

The constant operator 0a 

(a) is open at {a,b) iff a < pa; 

(b) it is open iff Pa > 0 A(pa) = oo. 

The simple operator 0a 
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(c) is open at {a,b) iffb< p^; 

(d) it is open iff Pa > 0 A(pa) = oo. 

Proof. ©A is open at (a, b) iff for some e > 0 we have A(a+£) < oo. This is clearly equivalent 
to a < Pa- 

Thus Pa > 0 and A(pa) < oo imply ©a is not open at (pa, b) for any 6 > 0, hence it is 
not open. Conversely if ©a is not open then Pa > 0 and A (a) < oo for some a,b > 0, but 
A(a + e) = cx) for any e > 0. This implies a = pa- 

The proof for the simple operator A{w) is similar. □ 


4.5. Operational closure of the set of open E. 


Lemma 37. Let a,b > 0. 

(a) The set of elementary operators open at (a, b) is closed under the arithmetical opera¬ 
tions of scalar multiplication, addition and multiplication. If E 2 is open at (a, b) and 
El is open at (a,E 2 (a, 6 )) then Ei(z, E 2 ( 2 :, tc)) is open at {a,b). 

(b) The set of open elementary operators is closed. 


Proof. Let c > 0 and let E, Ei, E2 be elementary operators open at {a,b). Then 
( 3 e > 0) E(a + e, 6 + e) < 00 ^ ( 3 e > 0) (cE)(a + e, 6 + e) <00 

( 3 ei > 0) El (a + ei, 6 + ei) < 00 and (3e2 > 0) E2 (a + 62, 6 + £2) < 00 

( 3 e > 0 )Ej (a + e, 6 + s) < 00 for z = 1 , 2 
(ds > 0 ) (Ei + E2) (o T £, 6 + e) < 00 
( 3 ei > 0) El (a + £1,6 + ei) < cxD and (^£2 > 0) E2 (a + £2, 6 + £2) < 00 

( 3 £ > 0 )Ej (o + £, 6 + £) < 00 for i = 1 , 2 
( 3 £ > 0) (E1E2) (a + £, 6 + £) < 00. 

Now suppose E 2 is open at (a, 6) and Ei is open at (a,Ei(a,6)). Then 

(3£2 > 0) E2(a + £2 ,5 + £2) < 00 and ( 3 £i > 0) Ei (a + £1, E2(o, 6) + £1) < 00 

^ 0 ) El (^a T £, E2(n T £^5 ^ T £) T £^^ ^ oo. 

This completes the proof for (a). Part (b) is proved similarly. 


□ 


The base operators that we will use as a starting point are the elementary operators 
E and all possible restrictions ©m of the standard operators © of combinatorics discussed 
below. More complex operators called composite operators will be fabricated from these base 
operators by using the familiar arithmetical operations of addition, multiplication, scalar 
multiplication and composition discussed in 11 14.21 


4.6. The standard operators on DOM[ 2 ;]. Following the lead of Flajolet and Sedgewick 
ra we adopt as our standard operators MSet (multiset). Cycle (undirected cycle), DCycle 
(directed cycle) and Seq (sequence), corresponding to the constructions by the same names.^ 

^Flajolet and Sedgewick also include Set as a standard operator, but we will not do so since, as mentioned 
in the second paragraph of §|21 for a given T, the series G(z, w) associated with Set(T) may very well not 
be elementary. For a discussion of mixed sign equations see §0 
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These operators have well known analytic expressions, for example, 

unlabelled multiset operator 1 + MSet(T) = exp 

labelled multiset operator MSet(T) = — 1 


4.7. Restrictions of standard operators. Let M C P. (We will always assume M is 
nonempty.) The 'M.-restriction of a standard construction A applied to a class of trees 
means that one only takes those forests in A(T) where the number of trees is in M. Thus 
MSet{ 2 , 3 }(^) takes all multisets of two or three trees from T. 

The Polya cycle index polynomials zi,..., Zm) are very convenient for expressing 
such operators; such a polynomial is connected with a permutation group H acting on an 
m-element set (see Harary and Palmer p. 35). For (T G H let aj be the number of 
j-cycles in a decomposition of a into disjoint cycles. Then 


Z (H, ,..., Zui ) 


1 

\H\ 


En 

<j£H j = l 


The only groups we consider are the following: 

(a) Sm is the symmetric group on m elements, 

(b) Dm the dihedral group of order 2m, 

(c) Cm the cyclic group of order m, and 

(d) Idm the one-element identity group on m elements. 

The M.-restrictions of the standard operators are each of the form Am := 
where A e {MSet, DCycle, Cycle, Seq} and Am is given by: 


operator unlabelled case 

MSetm(T) Z(Sm,T(z),...,T(z”^)) 
Cyclem(T) Z(Dm, T(z),..., T(z-)) 
DCyclem(T) Z(Cm, T(z),..., T(z™)) 
Seqm(T) Z(ldm,T(z),...,T(z™)) 


operator labelled case 

^^tm(T) (l/m!)T(z)™ 

C^em(T) {l/2m)T{zr 

[5^em(T) (l/m)T(z)- 
S;qm(T) T(zr 


Note that the labelled version of Am is just the hrst term of the cycle index polynomial 
for the unlabelled version, and the sequence operators are the same in both cases. We write 
simply MSet for MSetM if M is P, etc. 

In the labelled case the standard operators (with restrictions) are simple operators, whereas 
in the unlabelled case only Ap} and the SeqM ate simple. The other standard operators in 
the unlabelled case are not elementary because of the presence of terms T{z^) with j > 1 
when M ^ {!}. 


4.8. Examples of recursion equations. Table Q gives the recursion equations for the 
generating series of several well-known classes of trees. 


^m-flagged means one can attach any subset of m given flags to each vertex. This is just a colorful way 
of saying that the tree structures are augmented with m-unary predicates C/i,..., Um , and each can hold on 
any subset of a tree independently of where the others hold. 
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Recursion Equation 


Class of Rooted Trees 


w = z + zw 
w = z + 2 ;Seq(t(;) 
w = mz + mzSeq{w) 

7 /; 

w = ze 

w = z + 2 ;MSet(rt;) 
w = z + zMSet{2,3}(^i’) 
w = z + zSeq 2 iw) 
w = z + 2 ;MSet 2 (rc) 
w = z + zw‘^ 

w = z + z{w + MSet 2 (rc)) 
w = z + zMSetr(rc) 


chains 
planar 

m-flagged planar® 
labelled 
unlabelled 
unlabelled (0,2,3)- 
unlabelled binary planar 
unlabelled binary 
labelled binary 
unlabelled unary-binary 
unlabelled r-regular 


Table 1. Familiar examples of recursion equations 


4.9. Key properties of operators. Now we give a listing of the various properties of 
abstract operators that are needed to prove a universal law for recursion equations. The 
hrst question to be addressed is “Which properties does 0 need in order to guarantee that 
w = Q{w) has a solution?” 

4.10. Retro operators. There is a simple natural property of an operator 0 that guarantees 
an equation w = Q{w) has a unique solution that is determined by a recursive computation 
of the coefficients, namely 0 calculates, given T, the nth coefficient of 0(T) solely on the 
basis of the values of f(l),..., t{n — 1). 

Definition 38. An operator 0 is retro if there is a sequence a of functions such that for 
B = 0(A) one has bn = o'„(ai,..., a„_i), where Ui is a constant. 

There is a strong temptation to call such 0 recursion operators since they will be used to 
recursively dehne generating series. But without the context of a recursion equation there is 
nothing recursive about being a function of Oi,..., On-i ■ 

Lemma 39. A retro operator 0 has a unique fixpoint in DOM[^], that is, there is a unique 
power series T G DOM)^:] such thatT = 0(T). We can obtain T by an iterative application 
of 0 to the constant power series 0.' 

T = lim 0^(0). 

n—5-00 

7/0 is an integral retro operator then T G IDOM)^:]. 

Proof. Let a be the sequence of functions that witness the fact that 0 is retro. If T = 0(T) 
then 

t(l) = Ui 

t{n) = (T„(t(l),..., t(n — 1)) for n > 1. 

Thus there is at most one possible fixpoint T of 0; and these two equations show how to 
recursively find such a T. 

A simple argument shows that 0"^’''^(O) agrees with 0*^(0) on the first n coefficients, for 
all k > 0. Thus lim„^oo0”(O) is a fixpoint, and hence the hxpoint. If 0 is also integral 
then each stage 0”(O) G IDOM);^], so T G roOM)^:]. □ 
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Thus if 0 is a retro operator then the functional equation w = Q{w) has a unique solution 
T( 2 :). Although the end goal is to have an equation w = Q{w) with 0 a retro operator, for 
the intermediate stages it is often more desirable to work with weakly retro operators. 

Definition 40. An operator 0 is weakly retro if there is a sequenee a of funetions sueh that 
for B = 0(A) one has bn = , a„) . 


Lemma 41. 

(a) The set of retro operators is elosed. 

(b) The set of weakly retro operators is elosed and ineludes all elementary operators and 
all restrietions of standard operators. 

(c) If Q is a weakly retro operator then zQ and wQ are both retro operators. 

Proof. For (a), given retro operators 0, 0i, 02, a positive constant c and a power series T>0, 
we have 


M(ce)(T) 
[j“](ei + e2)(T) 

1^”] (SiS^XT) 

|j”| (Oioe^XT) 


c{iz"] e(T)) 

X“]e,(T) + X“]e2(T) 

n—1 

5 ;[X]ei(T)|*"-xeXT) 

i=i 

ff„(lX]e2(T),...,X”] e2(T)), 


where a is the sequence of functions that witness the fact that 0i is a retro operator. In 
each case it is clear that the value of the right side depends only on the hrst n — 1 coefficients 
of T. Thus the set of retro operators is closed. 

For (b) use the same proof as in (a), after changing the initial operators to weakly retro 
operators, to show that the set of weakly retro operators is closed. 

For an elementary operator E(; 2 , w) and power series T > 0 we have, after writing E(; 2 , w) 
as Ei>o 


X"]E(^,TW) = lz’']J2E,{z)T(zy 

j>0 

n 

j>0 i=0 

The last expression clearly depends only on the hrst n coefficients of T(; 2 ). Thus all elemen¬ 
tary operators are weakly retro operators. 

Let E(H, Zi,..., Zm) be a cycle index polynomial. Then for T G DOM);^] one has 

r nnrp/ _ j 0 if i does not divide n 

J “ \t(n/j) ifj|n. 

Thus the operator that maps T(; 2 ) to T{z^) is a weakly retro operator. The set of weakly 
retro operators is closed, so the operator mapping T to E(H, T(2:), ..., T(2;”^)) is weakly 
retro. Now every restriction Am of a standard operator is a (possibly inhnite) sum of such 
instances of cycle index polynomials; thus they are also weakly retro. 
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For (c) note that 


M(2e(T)) = |^"-■]e(T) 

n—1 

|2”](Te(T)) = ^«,|*”-i]e(T), 

i=l 

and in both cases the right side depends only on ti,..., tn-i- 

□ 


Lemma 42. 

(a) An elementary operator E(; 2 , ta) = retro ijf cqi = 0. 

(b) A restriction Am of a standard operator A is retro iff 1 

Proof. For (a) let T G DOM[; 2 ]. Then 

n 

j>0 i=0 

which does not depend on t(n) iff Coi = 0 . 

For (b) one only has to look at the dehnition of the Polya cycle index polynomials. □ 

The property of being retro for an elementary E(^, w) is very closely related to the neces¬ 
sary and sufficient conditions for an equation w = E{z, w) to give a recursive dehnition of a 
function T G DOM[; 2 ] that is not 0. To see this rewrite the equation in the form 

(1 ~ Goi)w = [ciqZ + 620 ^^ -|- ■ ■ ■) -|- (eiiz + ■ ■ -Jvj -h ■ ■ ■ . 

(We know that Coo = 0 as E is elementary.) So the hrst restriction needed on E is that 

Coi < 1. 

Suppose this condition on coi holds. Dividing through by 1 — coi gives an equivalent 
equation with no occurrence of the linear term z^w^ on the right hand side, thus leading to 
the use of coi = 0 rather than the apparently weaker condition coi < 1. 

To guarantee a nonzero solution we also need that Eo( 2 :) 7 ^ 0, and by the recursive con¬ 
struction these conditions suffice. 

Now that we have a condition, being retro, to guarantee that w = Q{w) is a recursion 
equation with a unique solution tc = T, the next goal is to hnd simple conditions on 0 that 
ensure this solution will have the desired asymptotics. 

4.11. Dominance between power series. It is useful to have a notation to indicate that 
the coefficients of one series dominate those of another. 

Definition 43. For power series A, B G DOM);^] we say B dominates A , written A < B , 
if Oj < bj for all j. 

Likewise for power series G, H G tc] we say H dominates G , written G < H , 

if dij < hij for alli,j. 

Lemma 44. The dominance relation < is a partial ordering on DOM);^] preserved by the 
arithmetical operations: for Ti, T 2 , T G DOM);^] and a constant c> 0, if Ti < T 2 then 

c ■ Ti < c ■ T 2 

Ti + T < T 2 + T 
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Proof. Straightforward. 


□ 


Ti T 

< 

T2 T 

Ti oT 

< 

T20T 

ToTi 

< 

T0T2 


4.12. The dominance relation on the set of operators. 

Definition 45. 

(a) For operators ©i, ©2 we say ©2 dominates ©1, symbolieally ©1 ^ ©2 , if for any T G 
D0M[2:] one has ©i(T) < ©2(T) . 

(b) For integral operators ©1, ©2 we say ©2 dominates ©1, symbolieally ©1 ©2 , if for 

any T G roOM[ 2 :] one has ©i(T) < © 2 (T) . 

As usual we continue our discussion mentioning only the general operators when the 
integral case is exactly parallel. It is straightforward to check that the dominance relation 
C is a partial ordering on the set of operators which is preserved by addition, multiplication 
and positive scalar multiplication. Composition on the right also preserves C, that is, for 
operators ©1 C ©2 and ©, 

©1 o © □ ©2 o ©. 

However composition on the left requires an additional property, monotonieity. 

The bivariate E in DOM[; 2 ,t(;] play a dual role, on the one hand simply as power series, 
and on the other as operators. Each has a notion of dominance, and they are related. 

Lemma 46. For E, F G DOM[ 2 :, tc] we have 

E < F ^ E □ F. 

Proof. Suppose E < F and let T G DOM[ 2 :]. Then 

ME(^,T) = 5^e.,|2”-]T(^r < Y.f„lz"-‘]T(zy = MF(2,T), 

SO 'E{z, T) < F( 2 :, T). As T was arbitrary, E □ F. □ 

4.13. Monotone operators. 

Definition 47. An operator © is monotone if it preserves <, that is, A < B implies 
©(A) < ©(B) for A, Be ©OM);^] . 

Lemma 48. //©i ©©2 and © is monotone then 

©o©i □ ©o©2. 


Proof. Straightforward. □ 

Lemma 49. The set of monotone operators is closed and includes all elementary operators 
and all restrictions of the standard operators. 


□ 


Proof. Straightforward. 
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4.14. Bounded series. 


Definition 50. For R > 0 let Aji{z) := J2n>i ^ series T G DOM[z] is hounded if 

T < Ar for some R > 0. 

An easy application of the Cauchy-Hadamard Theorem shows that T is bounded iff it is 
analytic at 0. 

The following basic facts about the series Ar{z) show that the collection of bounded series 
is closed under the arithmetical operations, a well known fact. Of more interest will be the 
application of this to the collection of bounded operators in Section 14.171 

Lemma 51. For c, R, Ri, R 2 > 0 


Ri F R2 


Ar^ < Ar.^ 

cAr 

< 

A(c+1)_R 

A_Ri + Ar^ 

< 

Ari+R2 

Ar^Ar.^ 

< 

Ari+R2 

Ar^ 0 Arj 

< 

A2(1+Ri+_R2)' 


Proof. The details are quite straightforward—we give the proofs for the last two items. 
(A„.Ah.)(^) = 

i>i i>i 


n>l i+j=n 

= E (E 

n>l i-\-j=n 

*.i>i 


< J2iRl + R2rz^ = ^Rr+R.{z)- 

n>l 


For composition, letting Rq = 1 + Ri + R 2 '. 


(Ahi o Ar2)(2;) 


(Aho 0 ArJ(z) 

= E 


(E^'-O’ 


i>l 


i>i 


')A 

E 

(E(-Ro^O^) 

*>i i>i 


i>l 

i>i 

n>l 

^2Ro^ ( 

z). 



□ 


4.15. Bounded operators. The main tool for showing that the solution tc = T to tc = 
Q{w) has a positive radius of convergence, which is essential to employing the methods of 
analysis, is to show that 0 is hounded. 
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Definition 52. For R> 0 define the simple operator Ar by 

A_r(w) = wfi 

i>i 

An operator 0 is bounded if {3R > O) (^Q{w) C Ar{z + w) j, that is, 

{3R > 0) (VT e roOM) (^0(T) < Ar{z + T)). 

Of course we will want to use integer values of R when working with integral operators. 

4.16. When is an elementary operator bonnded? The properties weakly retro and 
monotone investigated earlier hold for all elementary operators. This is certainly not the 
case with the bounded property. In this subsection we give a simple univariate test for being 
bounded. 

As mentioned before, any E G DOM[; 2 , w] plays a dual role in this paper, one as a bivariate 
power series and the other as an elementary operator. Each of these roles has its own 
dehnition as to what bounded means, namely: 

E<A^(2; + t(;) {'^i, j > l) (cij < [z" w^] Ar{z + w)'^ 

E T Ar{z + w) ^ (VT G D0M[;2]) T) < Ar{z + T)). 

The two dehnitions are equivalent. 

Lemma 53. Let E be an elementary operator. 

(a) Fi^z^w) is bounded as an operator iff Rii^z^z) is bounded as a power series. Indeed 

E( 2 :, tc) C Ar( 2 : + tc) ^ E{z, z) < A 2 r{z) forR>0 

^{z,z)<Ar{z) E{z,w)FAr{z+ w) forR>l. 

(b) The equivalence of bivariate bounded and operator bounded follows from 

R‘{z,w) < Ar{z + w) 'E{z,w) F Ar{z + w) forR >0 
E(;2, w) C Aij(2: + w) ^ 'E{z,w) < A2 r{z + w) forR>l. 

Proof. For (a) suppose R> 0 and E(; 2 , w) T Ar{z + w). Since 2 : G roOM[ 2 ;], we have 

Eiz,z) < = A2r(;^), 

i>i 

so Ri{z, z) is a bounded power series. 

Conversely, suppose R> 1 and Rii^z^z) < Aj:j(; 2 ). Then 

E(2,2) < 

i>i 

so for n > 1 

[R]E{z,z) < Rfi 

Then from E{z,w) = we have Cij < R^~^\ so 

E{z,w) < 

*j>i 
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<1 


*J >1 



Ar{z + w). 


Z^'W^ 


Applying Lemma UHl gives Fi{z, w) C Ar{z + w). 

For (b) the first claim is just Lemma EHl For the second claim suppose R > 1 and 
E(^, w) C Aji{z + w). From the hrst part of (a) we have £( 2 ;, z) < A2r{z) and then from the 
second part £( 2 ;, w) C A 2 r{z + w). 

□ 


Corollary 54. Given A G roOM[z], the constant operator ©a as well as the simple operator 
A{w) are bounded iff Pa > 0. 

4.17. Bounded operators form a closed set. 

Lemma 55. The set of bounded operators is closed. 


Proof. Let 0, ©i, ©2 be bounded operators as witnessed by the following: 

Q{w) ^Ar{z + w), Qi{w) C Aj:jj( 2 ; + «;) and © 2 (iy) ^AR.^i.^^'^)- With c > 0 we have from 
Lemma M 


(ce)(w) 

(©1 + ©2)(w) 

(0102) H 

(01 O ©2)(m;) 


□ cAr(z + w) O A(i+c)i?(^ + w) 

□ Ar^ (z + w) + Ar2 (z + w) O Ar^+r, (z + w) 

O Ar^(z + w)Ar2(z + w) o Ar^+r^(z + w) 

O Ar^(z + w) O Ar^(z + w) o A2 (i+r^+R2)2(z + w). 


□ 


Lemma 56. All restrictions of standard operators are bounded operators. 

Proof. Let A be a standard operator. Then for any M C P we have Am □ A, so it suffices 
to show the standard operators are bounded. But this is evident from the well known fact 
that 


MSet(t(;) □ Cycle(w) T DCycle(tc) 

T Sec\{w) = = Ai(tc) □ Ai( 2 : + tc). 

n>l 

So the choice of i? is i? = 1. □ 


4.18. When dominance of operators gives dominance of fixpoints. This is part of 
proving that the solution w = T to tc = 0(w) has a positive radius of convergence. 

Lemma 57. Let Tj satisfy the recursion equation T* = 0j(Tj) for i = 1,2 . If the 0j are 
retro operators, 0i T 02 , and 0i or ©2 is monotone then Ti < T 2 . 

Proof. Since each Qi{w) is a retro operator, by Lemma IHHl we have 

T, = lini ©,"(0). 

n—>cxD 

Let us use induction to show 

01^(0) < 02'^(O) 
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holds for n > 1. For n = 1 this follows from the assumption that 02 dominates ©i. So 
suppose it holds for n. Then 

0i”+^(O) < 0 i( 02 ''(O)) < 02 ”+^(O) if 01 is monotone 

0i”’''^(O) < 02 f 0 i"'(O)j < 02 "’'’^(O) if 02 is monotone. 

Thus Ti < T 2 . □ 

4.19. The nonzero radius lemma. To apply complex analysis methods to a solution T 
of a recursion equation we need T to be analytic at 0. 

Lemma 58. Let Q be a retro operator with 0(w) T + w). Then 

Q(w) C Aji(z + w)—Rw. 

Proof. Since 0 is retro there is a sequence of functions such that for T G DOM[^], 

[;^"]0(T) = an{tu... Pn-l). 

Let 

$(«;) ;= + 
n>2 

which is easily seen to be a retro operator. Choose such that for T G DOM[z] 

[z"]$(T) = a4ti,...,f„_i). 

Then, since Apti^z + w) = R{z + w) + ‘F(w), from the dominance of Q{w) by Ar{z + w) we 
have, for any ti > 0 and n> 2, 

(tl, . . . , ^ Rtji T , . . . , . 

As the left side does not depend on R we can put = 0 to deduce 

(fl ) • • • ) 1 ) ^ (^1) • • • ) l) y 

which gives the desired conclusion. □ 

Lemma 59. Let Q be a bounded retro operator. Then w = Q{w) has a unique solution 
w = T, and pt > 0. 

Proof. By Lemma EHl we know there is a unique solution T. Choose R > 1 such that 
Q{w) T Api{z + w). From Lemma IHHI we can change this to 

(22) Q{w) T Api{z + w) — Rw. 

The right side is a monotone retro operator, so Lemma EH says that the fixpoint S of 
Aji{z + w) — Rw dominates the hxpoint T of Q{w). Let 

S = Ar{z + S) - RS. 

To show Pt > 0 it suffices to show ps > 0. We would like to sum the geometric series 
Ar[z + S(z)); however since we do not yet know that S is analytic at z = 0 we perform an 
equivalent maneuver by multiplying both sides of equation (1^ by 1 — Rz — RS to obtain 
the quadratic equation 

{R + R^) + {R^z + Rz-l)S + Rz = 0. 
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The discriminant of this equation is 

D(^) = [R^z + Rz-lf - A{R + R^)Rz. 

Since D(0) = 1 is positive it follows that a/D(^) is analytic in a neighborhood of ^ = 0. 
Consequently 8 ( 2 ;) has a nonzero radius of convergence. 

□ 

4.20. The set of composite operators. The sets Oe and Oi of operators that we even¬ 
tually will exhibit as “guaranteed to give the universal law” will be subsets of the following 
composite operators. 

Definition 60. The composite operators are those obtained from the base operators, namely 

(a) the elementary operators E(z, w) and 

(b) the ML-restrictions of the standard operators: MSetM, Cycle^j, DCycle^ and Seq^j, 

using the variables z, w, scalar multiplication by positive reals, and the binary operations 
addition (- 1 -), multiplication (■) and composition (o). 

Lemma 61. The set of composite operators is closed under the arithmetical operations and 
all composite operators 0 are monotone and weakly retro. 

Proof. The closure property is immediate from the dehnition of the set of composite op¬ 
erators, the monotone property is from Lemma EHl and the weakly retro property is from 
Lemma El (b). □ 

An expression like z -|- 2 :Seq(r(;) that describes how a composite operator is constructed 
is called a term. Terms can be visualized as trees, for example the term just described and 
the term in dH) have the trees shown in Figure El (A small empty box in the hgure shows 
where the argument below the box is to be inserted.) Composite operators are, like their 





w 

Figure 5. Two examples of term trees 

counterparts called term functions in universal algebra and logic, valued for the fact that 
one has the possibility to ( 1 ) define functions on the class by induction on terms, and ( 2 ) 
one can prove facts about the class by induction on terms. 
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Perhaps the simplest explanation of why we like the composite operators 0 so much is: we 
have a routine procedure to convert the equation w = Q{w) into an equation w = E( 2 :,t(;) 
where E is elementary. This is the next topic. 

4.21. Representing a composite operator 0 at T. In order to apply analysis to the 
solution tc = T of a recursion equation w = Q{w) we want to put the equation into the form 
w = E(;2, w) with E analytic on T. The next dehnition describes a natural candidate for E 
in the case that 0 is composite. 

Definition 62. Given a base operator 0 and a T G ©OM[ 2 :] define an elementary operator 
E®’'^ as follows: 

(a) E^’T = E /or E an elementary operator. 

(b) For 0 = MSetM Ft E®’^ = ... ,T{z^)). 

(c) For0 = DCycleM /ef E®’^ = Y.^^^Z{Cm,w,T{z^),... ,T{zn). 

(d) For0 = CycleM /ef E®’^ = j:^^^Z{Dm,w,T{z^),... ,T{z^)). 

(e) For 0 = SeqM let E®’^ = 

Extend this to all composite operators using the obvious inductive definition: 

E^e,T _ 

E01+02,T ._ E®1’T _|_ E02,T 

E©i02,t _ e®1’Te®2-t 

E01O02,T ._ E®1’®2(T)^^ E02,T^ 

The dehnition is somewhat redundant as the Seq^^ operators are included in the elementary 
operators. 

Lemma 63. For 0 a composite operator and T G DOM[ 2 :] we have 

0(T) = E®’'^(;2,T). 

ITe will simply say that E®’”^ represents Q at T. 

Proof. By induction on terms. □ 

4.22. Defining linearity for composite operators. 

Definition 64. Let Q be a composite operator. We say 0 is linear (in w) if the elementary 
operator E®’^ representing Q at z is linear in w. Otherwise we say 0 is nonlinear (in w). 

Lemma 65. Let Q be a composite operator. Then the elementary operator E®’'^( 2 :,to) rep¬ 
resenting Q at T is either linear in w for all T G DOMM, or it is nonlinear in w for all 
T G D0M[;2]. 

Proof. Use induction on terms. □ 

4.23. When T belongs to DOM* [ 2 :]. 

Proposition 66. Let Q be a bounded nonlinear retro composite operator. Then there is 
a unigue solution w = T to w = Q{w), and T G DOM*[ 2 :], that is, px ^ (0, cxd) and 
T(pt) < 00. 
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Proof. From Lemma EHl we know that w = Q{w) has a unique solution T G DOM[z], 
and p '■= Pt > 0. Let E{z,w) be the elementary operator representing 0 at T. Then 
T = E(z, T). As 0 is nonlinear there is a positive coefficient of E with j > 2. Clearly 

T(a:) > eijX^T{xy for a: > 0. 

Divide through by T(x)^ and take the limsup of both sides as x approaches p~ to see that 
T(p) < cx), and thus p < oo. This shows T G roOM*[ 2 :]. □ 

4.24. Composite operators that are open for T. Many examples of elementary oper¬ 
ators enjoy the open property, but (restrictions of) the standard operators rarely do: only 
the various Seq^j and Ap} for A any of the standard operators. 

For the standard operators other than Seq, and hence for most of the composite operators, 
it is very important that we use the concept of ‘open at T’ when setting up for the Weierstrafi 
Preparation Theorem. 

Definition 67. Let T G DOM*[^]. A composite operator 0 is open for T is open 

at (p,T(p)). 

The next lemma determines when the base operators are open for a given T G DOM*[z]. 

Lemma 68. Suppose T G DOM*[ 2 ;] and let p G (0, oo) he its radius of convergence. Then 
the following hold: 

(a) An elementary operator E is open for T iff it is open at (p, T(p)). 

(b) A constant operator Qa{w) is open for T iff p < p^. 

(c) A simple operator A{w) is open for T iff T(p) < pA- 

(d) Seqj^ is open for T iff M is finite or T{p) < 1. 

(e) MSetM is open for T iff M = {1} or p < 1. 

(f) DCycleM, or Cycle^, is open for T iff M = {1} or (M is finite and p < 1) or (M is 
infinite and p, T(p) < 1). 

Proof. For (a) note that an open operator represents itself at T. For (b) and (c) use Lemma 
IMl For (d) note that Seqf^(t(;) is the simple operator A{w) := (‘^) applies. 

For (e) let E := E®’”^ where 0 := MSetM- Then 

E{z,w) := 

mGM 

If M = {1} then E{z, w) = w and (c) applies. So suppose M y {!}. The term T{z^) appears 
in E{z, w), and this diverges at p -f e if p > 1. Thus p < 1 is a necessary condition for E to 
be open for T. 

So suppose p < 1. The representative for MSet dominates the representative of any MSetM- 
Thus for any x G (0, ffp) and p > 0: 

E(x, y) < exp (S < cx). 

m>2 

Since one can hnd £ > 0 such that the right hand side is hnite at {p + e, T(p) -1 -e), it follows 
that MSetM is open for T when p < 1. 
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For (f) let E := E®’’^ where 0 := DCyclej^. Then 

E(z,u.) := J]z(C„.k,,T(22)....,T(2”)) 

m€M 



A{w) B(z) 


If M = {1} then, as before, there are no further restrictions needed as E(^, w) := w. So now 
suppose M 7 ^ {!}. The presence of some T{z^) with fc > 2 in the expression for E(z,t(;) 
shows, as in (e), that a necessary condition is p < 1. This condition implies pe > ^/p■ 

If M is hnite then pa = oo, and pe > y^, consequently E is open at (p, T(p)). 

If M is inhnite then pA = 1. Suppose E is open at (p, T(p)). Then A(T(p) +e) converges 
for some £ > 0, so T(p) < 1. The conditions p, T(p) < 1 are easily seen to be sufficient in 
this case. 

For the Cyclcj^ case let E := E®’’^ where 0 := Cycle^. 


Cycles (T(^)) 


Thus 


E{z, w) 


^DCyc\e^{T{z)) 

is odd 
is even. 


4 ^ 


2T(;2)T(;22)(m-l)/2 

T(;2)2T(^2)("^-2)/2 + T{z^)'^/‘^ iim 


mSM 


1 

2 



mEM 



v?(fc) 

k 




2 wT (^ 2 )(»^- 2)/2 

y^2T(^2)(m-l)/2 


+ T{z^)^/‘^ 


if m is odd 
if m is even 


and we can use the same arguments as for DCycle. 


□ 


4.25. Closure of the composite operators that are open for T G DOM*[z]. 

Lemma 69. Suppose T G DOM*[ 2 :]. Then the following hold: 

(a) The set of composite operators that are open for T is closed under addition, scalar 
multiplication and multiplication. 

(b) Given composite operators 0i, 02 with 02 open for T and 0i open for Ti := 02(T), 
the composition 0i o 02 is open for T. 

Proof. Just apply Lemma EH □ 


4.26. Closure of the composite integral operators that are open for T G IDOM*[z]. 

Definition 70. roOM*[z] = IDOM[z] fl DOM* [z]. 

Lemma 71. Suppose T G IDOM*[z]. Then the following hold: 
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(a) The set of integral composite operators that are open for T is closed under addition, 
positive integer scalar multiplication and multiplication. 

(b) Given integral composite operators ©i, ©2 with ©2 open for T and ©1 open for Ti := 
© 2 (T), the composition ©1 o ©2 is integral and open for T. 

Proof. This is just a repeat of the previous proof, noting that at each stage we are dealing 
with integral operators acting on IDOM[z]. □ 

4.27. A special set of operators called O. This is the penultimate step in describing the 
promised collection of recursion equations. 

Definition 72. Let O he the set of operators that can he constructed from 

(a) the hounded and open elementary operators E(z, tc) and 

(b) theMi-restrictions of the standard operators: MSetM, Cycles, DCyclej^ andSeq-^, where 
in the case of the cycle constructions we reguire the set M to he either finite or to satisfy 
EmGMilM = 00. 

using the variables z, w, scalar multiplication hy positive reals, and the binary operations 
addition (+), multiplication (■) and composition (o). 

Within O let Oe be the set of bounded and open elementary operators; and let Oj be 
the closure under the arithmetical operations of the hounded and open integral elementary 
operators along with the standard operators listed in (h). 

Clearly (9 is a subset of the composite operators. 

Lemma 73. 

(a) Every Q E O is a bounded monotone and weakly retro operator. 

(b) Each of the sets O, Oe, Oj is closed under the arithmetical operations. 

Proof. For (a) we know from our assumption on the elementary operators in O and Lemma 
l^that the base operators in O are bounded—then Lemma shows that all members of O 
are bounded. All members of O are monotone and weakly retro by Lemma IHTl Regarding 
(b), use Lemma IT7I fbl for Oe, and Definition 1721 for the other two sets. □ 

Lemma 74. Let & E Oj. // T G roOM*[ 2 ;] and ©(T)(pt) < cxo then © is open for T. 

Proof. Since T E IDOM* we must have p := Pt < 1- Let 

O* := {OeOj: ©(T)(p) < 00}. 

An induction proof will show that for Q E O* we have © open for T. The elementary base 
operators of O* are given to be open, hence they are open for T. The restrictions of the 
standard operators in O* are covered by parts (d)-(f) of Lemma IHHl with one exception. 
We need to verify in certain DCycle and Cycle cases that T(/9) < 1. In these cases one has 
M infinite, and then we must have T(p) < 1 in order for ©(T) to converge a.t z = p since 
EmewilM = 00. 

For the induction step simply apply Lemma ITTl □ 

38 


4.28. The Main Theorem. The following is our main theorem, exhibiting many 0 for 
which w = Q{w) is a recursion equation whose solution satisfies the universal law. Several 
examples follow the proof. 


Theorem 75. Let 0i be a nonlinear retro member of Oe, respectively Oi, and let A(z) G 
DOM[z], respectively A(z) G roOM[z], be such that A(pa) = oo. Then there is a unique 
T G D0M[2;], respectively T G IDOM[z], such that T = A( 2 :) + 0i(T). The coefficients of 
T satisfy the universal law (*) in the form 


t{n) 


pE,{p,T{p)) 
27rE^^(p, T(p)) 


■ p 


for n = d mod q. 


Otherwise t{n) = 0. Thus (*) holds on {n : t{n) > 0}. The constants d, q are from the shift 
periodic form T{z) = z^Y[z^). 


Proof. Let Q{w) = A( 2 :) + Qi{w), by Lemma [721 a member of Oe, respectively Oj. By 
Proposition IHHl there is a unique solution w = T to tc = A(z) + 0(w) and T G DOM*[ 2 ;]. 
Let = E®’’^. Then the elementary representative E of A(^) + Q{w) is given by 

E{z.,w) := A(2:) + Ei(z,rc). 

We will verify the hypotheses (a)-(e) of Theorem | 2 H 1 

T = E(^, T) by LemmaESl this is| 2 Hl(a). The fact that T G DOM* [2;] is| 2 Hl(b). By Lemma 
inSlwe get|^(c). Since A(0) = 0 and A 7^ 0 it follows that A^ 7^ 0. As E(z, 0) = A{z) it 
follows that E 2 7^ 0. This isEHl fdL 

To show 0 is open for T we note that in the case of the operators coming from Oe they 
are given to be open elementary operators; and for the case they are coming from Oj use 
Lemma ITil This gives | 2 H 1 (e) • □ 


4.29. Applications of the main theorem. One readily checks that all the recursion equa¬ 
tions given in Table Q satisfy the hypotheses of Theorem [721 One can easily produce more 
complicated examples such as 

w = 3z^ + 2 ;^Cycle( w) + w ^DCycle(w) -f MSet 2 (tc). 

Such simple cases barely scratch the surface of the possible applications of Theorem [721 
Let us turn to the more dramatic example given early in iD, namely: 

w = z + zMSet(^Seq( ^ ^ (2" + 1)(DCyclep,i,^^ 3 (M;))''. 

nGOdd nGEven 


We will analyze this from ‘the inside out’, naming the operators encountered as we work up 
the term tree. First we give names to the nodes of the term tree: 


$4 := MSet($3) <F5 

A(^) := 2; 01 


DCyclep.i„„(to) $3 — Seq(4i) 
E«Ev.„( 2” + I)*"" I’s := <1>5(4>4) 

2:<F4<F6 


Now we argue that each of these operators is in Oj: 

(a) $1 is an elementary (actually simple) integral operator with radius of convergence 1/6. 
Thus it is bounded. Since it diverges at its radius of convergence, it is open. Thus 
$1 G Oj. 
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(b) $2 is a restriction of DCycle to the set of prime numbers; since XlmePrimes 1/?^ = oo 
we have $2 G C)/- 

(c) $3 is in Oi as it is a composition of two operators in Oi. 

(d) $4 is in Oj as it is a composition of two operators in Oi. 

(e) $5 is an elementary (actually simple) integral operator with radius of convergence 1 / 2 . 
Thus it is bounded. Since it diverges at its radius of convergence, it is open. Thus 
$5 G Oj. 

(f) $6 is in Oi as it is a composition of two operators in Oi. 

(g) 01 is in Oi as it is a product of two operators in Oi. 

(h) 01 is a nonlinear retro operator in Oi. 

Thus we have an equation w = A(^) + 0i(w) that satishes the hypotheses of Theorem 
17^ consequently the solution w = T{z) has coefficients satisfying the universal law. 

4.30. Recursion specifications for planar trees. When working with either labelled 
trees or planar trees the recursion equations are elementary. Here is a popular example that 
we will examine in detail. 

Example 76 (Planar Binary Trees). The dehning equation is 


w = z zw'^. 


This simple equation can be handled directly since it is a quadratic, giving the solution 


Clearly p = 1/2 and for n > 1 we have t{2n) = 0, and Lemma ^1 gives 



For illustrative purposes let us examine this in light of the results in this paper. Note that 

E(2:,tc) := z + zw"^ 

is in the desired form A( 2 :) + 0i(tc) with A(pa) = 00 and 0i a bounded retro nonlinear 
(elementary) operator. 

The constants d, q of the shift periodic form are given by: 

(a) d = 1 as Eo( 2 :) = 2 : implies Eq = {1}. 

(b) g = 2 as Eq = {1}, E 2 = {1}, and otherwise Ej = 0; thus 1JE„ + (n — l)d = 
{Eq - 1) U (E 2 + 1) = {0, 2}, so q = gcd{0, 2} = 2. 

Thus t{n) > 0 implies n = 1 mod 2, that is, n is an odd number. For the constant in the 
asymptotics we have 


^^{z, w) = 1 + 
E'ujw{zj 2z. 


In this case we know p = 1/2 and T(p) = 1 (from solving the quadratic equation), so 


E,(p,T(p)) = 2 
E^^„(p, T(p)) = 1. 
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Thus 


t{n) 


pEz(p, T(p)) ^ 


27rEu,t„(p, T(p)) 

V 27r 

— ■ 2"'n“^/^ for n = 1 mod 2. 

TT 


4.31. On the need for integral operators. Since the standard operators, and their re¬ 
strictions, are dehned on DOM[z] it would be most welcome if one could unify the treatment 
so that the main theorem was simply a theorem about operators on DOM[z] instead of hav¬ 
ing one part for elementary operators on DOM[z], and another part for integral operators 
acting on roOM[2;]. However the following example indicates that one has to exercise some 
caution when working with standard operators that mention T(x-^) for some j > 2. 

Let 

Q{w) := |(l-F MSet2(tc)). 

This is 1/2 the operator one uses to dehne (0,2)-trees. This operator is clearly in O and 
of the form A(z) -I- Qi{w)] however it is not in either Oe or Oj, as required by the main 
theorem. 

0 is clearly retro and monotone. Usual arguments show that w = Q{w) has a unique 
solution w = T which is in DOM*, and we have 

(23) T(p) = Ip + + ipT(p2). 

Since 0(T) involves T(z^) it follows that p < 1 (for otherwise T(p^) diverges). 

Suppose p < 1. Following Polya let us write the equation for T as w = E(2:, tc) where 

E{z,w) := ^z + ^zw^ + ^zT{z^). 

Then the usual condition for the singularity p is 1 = E^(p, T(p)), that is 

(24) 1 = i(p2T(p)) = 

SO pT(p) = 2. 

Putting T(p) = 2/p into equation (1^ gives 

2 1 1 1 9x 

“ — TjP -^ 

P 2 p 4 

so 

4 = 2p2 + p2T(p2). 

Since T(p^) < T(p) = 2/p we have 

4 < 2p^ + 2p, 


a contradiction as p < 1. 

Thus p = 1, and we cannot apply the method of Polya since E{z,w) is not holomorphic 
at (1,T(1)). 
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5. Algorithmic Aspects 


5.1. An algorithm for nonlinear. Given a term that describes a composite oper¬ 

ator 0 there is a simple algorithm to determine if 0 is nonlinear. Let us use the abbreviation 
A for the various standard unary operators and their restrictions as well as the elementary 
operators E(; 2 , w). We can assume that any occurrence of a M-restriction of a standard op¬ 
erator A in $ is such that M ^ {1} since if M = {1} then Am is just the identity operator. 
Given an occurrence of a A in $ let Ta be the full subtree of <h rooted at the occurrence of 
A. 

An algorithm to determine if a composite 0 is nonlinear 

• First we can assume that constant operators 0 a are only located at the leaves of the 
tree. 

• If there exists a A in the tree of <h such that a leaf w is below A, where A is either a 
restriction of a standard operator or a nonlinear elementary B, then 0 is nonlinear. 

• If there exists a node labelled with multiplication in the tree such that each of the 
two branching nodes have a tc on or below them then 0 is nonlinear. 

• Otherwise 0 is linear in w. 

Proof of the correctness of the Algorithm. (A routine induction argument on terms.) □ 

6. Equations w = G[z,w) with mixed sign coefficients 

6.1. Problems with mixed sign coefficients. We would like to include the possibility of 
mixed sign coefficients in a recursion equation w = G{z,w). The following table shows the 
key steps we used to prove (*•) holds in the nonnegative case, and the situation if we try the 
same steps in the mixed sign case. 



G G R^'J[[; 2 ,w]] 
Nonnegative G 

G G K[[z,w]] 

Mixed Signs G 

Property 

Reason 

Reason 

(3!T) {T = G{z, T) 

p > 0 

p < 00 

T(p) < cx) 

G holomorphic in nbhd of T 
G^^,^(p, T(p)) 7 ^ 0 

G^p.TW) 7^0 

DomSing = {z : z'^ = p‘^} 

9oi = 0 

G is bounded 

G is nonlinear in w 

G is nonlinear in w 

G(p -(- £, T(p) -|- e:) < CX) 
G is nonlinear in w 

Go( 2 :) 7 ^ 0 

Spec Gy, ( 2 ;, T(; 2 )) is nice 

9oi = 0 

G is abs. bounded 
(?) 

(?) 

(?) 

(?) 

(?) 

(?) 


As indicated in this table, many of the techniques that we used for the case of a nonnegative 
equation do not carry over to the mixed case. 

(a) To show that a unique solution tc = T exists in the mixed sign case we can use the 
retro property, precisely as with the nonnegative case. The condition for G G tc]] 
to be retro is that 5^01 = 0. 

(b) To show p > 0 in the nonnegative case we used the existence of an i? > 0 such 
that E(2:, T) < Aji(^z -|- T). In the mixed sign case we could require that G{z, T) be 
absolutely dominated by Aji[z + T). 
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(c) To show p < oo and T(p) < cx) we used the nonlinearity of E(2:,ta) in w. Then 
T = E(z,T) implies T(x) > CijX^T^xY for some e^- > 0 with j > 2. This conclusion 
does not follow in the mixed sign case. 

(d) After proving that T G DOM*[z], to be able to invoke the theoretical machinery of §01 
we required that E be open at (p, T(p)), that is, 

(ds > O) T(p) + £■) < cx)^. 

This shows E is holomorphic on a neighborhood of T. In the mixed sign case there 
seems to be no such easy condition unless we know that J2ij 

(e) In the nonnegative case, if E is nonlinear in w then E^^ does not vanish, and hence 
it cannot be 0 when evaluated at (p, T(p)). With the mixed signs case, proving 
Gu,i„(p, T(p)) 7^ 0 requires a fresh analysis. 

(f) A similar discussion applies to showing G2(p, T(p)) ^ 0. 

(g) Finally there is the issue of locating the dominant singularities. The one condition we 
have to work with is that the dominant singularities must satisfy Gw{z, T(^)) = 1. In 
the nonnegative case we were able to use the analysis of the spectrum of E^: 

Spec(^E^( 2 :,T( 2 ;)) j = (J + (n - 1) 0 T. 

n 

This tied in with an expression for the spectrum of T(z)). However for the mixed 
case we only have 

Spec(^E^(2;, T(2;)) j C En + {n - 1 ) Q T. 

n 

In certain mixed sign equations one has a promising property, namely 

Gw{z,T{z)) G DOM[^]. 

This happens with the equation for identity trees. In such a case put G^ [z, T(2;)) in its 
pure periodic form U(z^). Then the necessary condition on the dominant singularities 
becomes simply 

= pP. 

If one can prove p = q, as we did with elementary recursions, then the dominant 
singularities are as simple as one could hope for. 

There is clearly considerable work to be done to develop a theory of solutions to mixed sign 
recursion equations. 

6.2. The operator Set. The above considerations led us to omit the popular Set operator 
from our list of standard combinatorial operators. In the equation w = z + zSet{w) for 
the class of identity trees (that is, trees with only one automorphism), one can readily show 
that the only dominant singularity of the solution T is p. But if we look at more complex 
equations, like 

w = z -\- z^ -\- z^ -\- zSet(Set(tc)MSet(tc)), 

the difficulties of determining the locations of the dominant singularities appear substantial. 
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Example 77. Consider the restrictions of the Set operator 
SetM(T) = ^ Setm(T), where 

mgM 

Set„(T) = 

Thus in particular 

Set2(T) = l(T{zf-T:{z^)). 

The recursion equation 

w = z + Set2(tc) 

exhibits different behavior than what has been seen so far since the solution is T(;2) = 
which is not a proper inhnite series. The solution certainly does not have coefficients 
satisfying the universal law, nor does it have a hnite radius of convergence that played such 
an important role. 

We can modify this equation slightly to obtain a more interesting solution, namely let 

Q{w) = z + z^ + Set2(tc). 

Then 0 is integral retro, and the unique solution tc = T to tc = 0(w) is 

T(z) = z + z^ + z^ + z^ + 2z^ + 3z^ + 6z'^ + 11;^® + ■ ■ ■ 
with t(n) > 1 for n > 1. Consequently we have the radius of convergence 
(25) p:=p^e[0,l]. 


We will give a detailed proof that T has coefficients satisfying the universal law, to hint at 
the added difficulties that might occur in trying to add Set to our standard operators. 

Let 

0i(tc) = z + z'^ +-w'^ 

a bounded open nonlinear retro elementary operator, hence an operator in Oe to which 
the Main Theorem applies. For A G TOOM note that 0(A) < 0i(A), so we can use the 
monotonicity of 0i to argue that 0'^(O) < 0i”(O) for all n > 1. Thus T is dominated by 
the solution S to tc = Qi{w). At this point we know that Pt > Ps > 0. 

Since p G (0,1] we have T{x‘^) < T(a;) for x G (0,p). Thus for x G (0,p) 

T(a;) > X + x'^+ ^[T{x)'^-T{ x)), 


or 

^T(a;) > X + x'^ + ^T{x)'^. 

Thus T(a;) cannot approach cx) as a; —> p“. Consequently T(p) < oo, and then we must also 
have p < 1. By dehning 


G{z, w) 


z + z^ + ^ (w^ 
1 2 1 ■ 






we have the recursion equation w = G(z,w) satished by tc = T, and G(z,w) has mixed 
signs of coefficients. 
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As p < 1 we know that G{z,w) is holomorphic in a neighborhood of the graph of T, so 
a necessary condition for 2 ; to be a dominant singularity is that Gw{z,T{z)) = 1, that is 
T(z) = 1. Since T is aperiodic, this tells us we have a unique dominant singularity, namely 
z = p, and we have T(p) = 1. 

Differentiating the equation 

T(z) = ^ - tT(z^} 

gives 

T{z) = l + 2z + T{z)T\z) - zT'{z‘^) 

or equivalently 

(1 - T(; 2 ))T'(; 2 ) = (1 + 2z) - zT'{z‘^) for |^| < p. 

Since p < 1 we know that 

Mm ((1 + 2 j) - ^T'(j=)) = (l + 2p)-pTV)^ 

|z|<p 

Let A be this limiting value. Consequently 

(26) lim (1 - T(j))T'( 2 ) = A. 

\z\<p 

By considering the limit along the real axis, as a: —> p“, we see that A > 0, so 

(1 + 2p) - pr(p2) = A > 0. 

Let 

F(2;,w) := w - (^z + z'^ + - ^T{z'^)^. 

Then 

F,{z,w) = -(l + 2z- zT\z‘^)^ = zT\z^)-{l + 2z), 
so 

F4p,T(p)) = pTV)-(l + 2p) = -A. 

If A > 0 then F,j(p, T(p)) < 0; and since F^^ = —2 we have 

f4p,t(p))f^^(p,t(p)) > 0. 

This means we have all the hypotheses needed to apply Proposition [TT] to get the square 
root asymptotics which lead to the universal law for T. 

To conclude that we indeed have the universal law we will show that A > 0. Let a G [p, 1]. 
Then for x G (0,p) we have T(a;^) < Q;T(a:), and thus for x G (0,p) 

T(a;) > a; + ^ ^T(a;)^ — aT(a;)j. 

Let 

U(a:) = X + a:^ + ^ ^U(x)^ — aU(x)j. 

U(x) = ^ ^(2 + a) - a /(2 + aY - 8(x + x^) j 
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Then 



Pu = “2 + 4 V4T^(2T"^. 

Now for X G / ;= (0, min(p, pu)) 

(2 + a)T(x) - T(x) 2 > (2 + a)U(x)-U(x)2 
so 

U(x)2-T(x) 2 > (2 + a)(U(x)-T(x)). 

Thus U(x) ^ T(x) for x G /. If U(x) > T(x) on I then 

U(x) + T(x) >2 + 0 ! for X G/. 

But this is impossible since on I we have 

U(x) < U(pu) = l + a/2 
T(x) < T(p) = 1. 

Thus we have 

U(x) < T(x) on I. 

If Pu < P then U(p; 7 ) < T(p[/), which is also impossible. Thus 

P < Pu- 

Now dehne a function / on [p, 1] that maps o G [p, 1] to pu as given in the preceding lines, 
that is: 

/(o) = (—1/2) + (l/4)^4 + 2(2 + o)^. 

Then o G [p, 1] implies /(o) G (p, 1]. Calculation gives /^(l) = 0.536..., so 

p < 0.54 

Since ps = (v^ — l)/2 = 0.366 .. . we have p^ < ps, and then 

r(p2) < s'{p% 
so 

-A = pTV)-(l + 2p) < pSV)-(l + 2p) < 0 
since xS'(x^) — (1 + 2x) < 0 for x G (0, 0.55). This proves A > 0, and hence the universal 
law holds for the coefficients of T. 

This example shows that the generating function T* for the class of identity (0,l,2)-trees 
satishes the universal law. We have T* dehned by the equation w = z + z * Set{i^ 2 }('?i'), and 
it turns out that = fn+i- (We discovered this connection with T* when looking for the 
hrst few coefficients of T in the On-Line Encyclopedia of Integer Sequences.) 

Example 78. In the 20 Steps paper of Harary, Robinson and Schwenk HZI the asymp¬ 
totics for the class of identity trees (those with no nontrivial automorphism) was successfully 
analyzed by hrst showing that the associated recursion equation 

w = z zSet{w) 

has a unique solution re = T G roOM*[ 2 :]. Then 

G{z,w) := 2 -exp (^(-l)™+^T(^’")/m) 

m>2 
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is holomorphic in a neighborhood of the graph of T. One has 

z + Gw{z,w) = G{z,w), 

so the necessary condition Gw{z, T(; 2 )) = 1 for a dominant singularity is just the condition 
T{z) = 1 + z. p is the only solution of this equation on the circle of convergence as T — > 0 

and is aperiodic. Consequently the only dominant singularity is z = p. 

The equation for identity trees is of mixed signs. From 

(27) Tiz) = z'^{\ + z>f 

i>i 

we can calculate the first few values of t{n) for identity trees:® 


t{l) 

«(2) 

*(3) 

t(4) 

((5) 

t(6) 

t{7) 

((8) 

m 

((10) 

1 

1 

1 

2 

3 

6 

12 

25 
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Returning to the dehnition of G we have 

G{z,w) = z + ze^^(^^{-ir+\t{l)z^+ t{2)z^^+ ■■■)/mY jnl 

n>0 m>2 

= z + ze^^(^ - {z^ + z"^ + ■ ■ ■)/2 + {z^ + z^ + ■■■)/?, 

n>0 

- [z^^ + z^ + ■■■)/jnl 
= z + ze^'^(^- Z^I2 + - 3^74 + ■ ■ ■ )”/n! 

n>0 

= (^1 - ;2V2 + ^^73 - 5//8 + ■ ■ ■ ) • 

Thus for some of the z^w^ the coefficients are positive, and some are negative; G{z,w) is a 
mixed sign operator. 

If one were to form more complex operators 0 by adding the operator Set and its re¬ 
strictions SetM to our set of Standard Operators, then there is some hope for proving that 
one always has the universal law holding for the solution to tc = 0(tc), provided one has a 
solution that is not a polynomial. 

The hope stems from the fact that although the G{z,w) associated with 0 may have 
mixed sign coefficients, when it comes to the condition T(; 2 )) = 1 on the dominant 

singularities we have the good fortune that G^(; 2 , T( 2 :)) > 0, that is, it expands into a 
series with nonnegative coefficients. The reason is quite simple, namely using the bivariate 
generating function we have 

Set„(T) = [m™] exp (^^(-l)"-VT(a;’")/n). 

n>l 

Letting Gm{z,w) be Zm(Sm, if, —T( 2 :^),..., (—we have 

Gmiz:w) = K]e“"'-exp(5^(-ir-VT(a;")/n), 

n>2 

®One can also look up sequence number A004111 in the On-Line Encyclopedia of Integer Sequences. 
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thus 


u . _ s 

= g-K]^“"'exp(5:(-ir‘«”T(x”)/„) 

n>2 

= K] ■ exp 

n>2 

= ■ exp ^(-l)”-iM^T(x’^)/n) 

n>2 

= Gm-liz,w). 

Consequently if we put G^( 2 ;,T(^)) into its pure periodic form \J{z'p) then we have the 
necessary condition z'^ = fP on the dominant singularities z. Letting z^\{z'^) be the shift 
periodic form of T(z) it follows that q\p. If we can show that p = q then DomSing = {z : 
z^ = p'^}, which is as simple as possible. Indeed this has been the case with the few examples 
we have worked out by hand. 

7. Comments on Background Literature 

Two important sources offer global views on hnding asymptotics. 

7.1. The “20 Step algorithm” of [T7j. This 1975 paper by Harary, Robinson and Schwenk 
is in good part a heuristic for how to apply Polya’s method^°, and in places the explanations 
show an affinity for operators close to the original ones studied by Polya. For example it 
says that G{z,w) should be analytic for \z\ < ^/pr and |tc| < oo. This strong condition on 
w fails for most of the simple classes studied by Meir and Moon, and hence for the setting 
of this paper. 

The algorithm of 20 Steps also discusses how to find asymptotics for the class of free trees 
obtained from a rooted class defined by recursion. Given a class T of rooted trees let U be 
the associated class of free (unrooted) trees, that is, the members of U are the same as the 
members of T except that the designation of an element as the root has been removed. Let 
the corresponding generating series^^ be T(z) and U(z). 

The initial assumptions are only two: that T is not a polynomial and it is aperiodic. Step 
2 of the 20 steps is: express U( 2 :) in terms of T( 2 ;) and T(z^). 20 Steps says that Otter’s 
dissimilarity characteristic can usually be applied to achieve this. Step 20 is to deduce that 
u{n) ~ Cp~’^rr^l‘^. 

This outline suggests that it is widely possible to hnd the asymptotics of the coefficients 
u{n), and evidently this gives a second universal law involving the exponent —5/2 instead 
of the —3/2. Our investigations suggest that determining the growth rate of the associated 
classes of free trees will be quite challenging. 

Suppose T is a class of rooted trees for which the Polya style analysis has been successful, 
that we have found the radius of convergence p G (0,1), that T(p) < oo, and that t{n) ~ 

^'^The paper also has a proof that the generating function for the class of identity trees (defined by 
w = z + zSet(iy)) satisfies the universal law. 

^^20 Steps uses t(n) to denote the number of rooted trees in T on n + 1 points, whereas u{n) denotes the 
number of free trees inhl on n points. We will let t[n) denote the number of free trees in T on n points, as 
we have done before, since this will have no material effect on the efficacy of the 20 Steps. 
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C. What can we say about the generating function U(^) for the corresponding class 
of free trees? 

Since we have the inequality t{n)/n < u{n) < t{n) (note that one has only n ways to 
choose the root in a free tree of size n), it follows by the Cauchy-Hadamard Theorem that 
U(z) has the same radius of convergence p as T{z). From this and the asymptotics for t{n) 
it also follows that one can hnd Ci, (72 > 0 such that for n > 1 

Cip-^n-^/^ < u{n) < C 2 p-"n-^^^. 

Thus u{n) is sandwiched between a —5/2 expression and a —3/2 expression. In the case that 
T is the class of all rooted trees, Otter 1211 showed that 

UW = ^(TW-Set2(T)) = ^(tW-1(TW"-T(2"))), 

and from this he was able to hnd the asymptotics for u{n) with a —5/2 exponent. 

However let T be the class of rooted trees such that every node has either 2 or 5 descending 
branches. The recursion equation for T(z) is 

T(z) = z + z(MSet2(T) + MSet5(T)), 

and T(z) is aperiodic. By Theorem I7H1 we know that the coefficients of T(^) satisfy the 
universal law 

t{n) - Cp-^n-^^^. 

Let U be the corresponding set of free trees. Note that when one converts a rooted tree 
T in T to a free tree F, a root with 2 descending branches will give a node of degree 2 in 
F, and a root with 5 descending branches will give a node of degree 5 in F. Any non root 
node with 2 descending branches will give a node of degree 3 in F; and any non root node 
with 5 descending branches will give a node of degree 6 in F. Thus F will have exactly one 
node of degree 2 or degree 3, and not both, so one can identify the node that corresponds 
to the root of T. This means that there is a bijection between the rooted trees on n vertices 
in T and the free trees on n vertices in U. Consequently t{n) = u{n), and thus 

u{n) ~ 

Clearly u{n) cannot also satisfy a —5/2 law. Such examples are easy to produce. 

Thus it is not clear to what extent the program of 20 Steps can be carried through for 
free trees. It seems that free trees are rarely dehned by a single recursion equation, and 
it is doubtful if there is always a recursive relationship between U( 2 :) and T(z), T( 2 :^),.... 
Furthermore it is not clear what the possible asymptotics for the u{n) could look like—is 
it possible that one will always have either a —3/2 or a —5/2 law? Since a class U derived 
from a T which has a nice recursive specihcation can be dehned by a monadic second order 
sentence, there is hope that the u{n) will obey a reasonable asymptotic law. (See Q5 in §|H1) 
In 20 Steps consideration is also given to techniques for calculating the radius of conver¬ 
gence p of T and the constant C that appears in the asymptotic formula for the t{n). In 
this regard the reader should consult the paper of Plotkin and Rosenthal I2nj as there are 
evidently some numerical errors in the constants calculated in 20 Steps. 
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7.2. Meir and Moon’s global approach. In 1978 Meir and Moon [20] considered classes 
T of trees with generating functions T(z) = 

(1) t(l) = 1; 

(2) T can be obtained by taking certain forests of trees from T and adding 
a root to each one (this choice of certain forests is evidently a ‘construc¬ 
tion’); 

(3) this ‘construction’ and ‘conditions implicit in the dehnition’ of T give 
rise to a ‘recurrence relation’ for the t{n), evidently a sequence a of 
functions such that t{n) = cT„(t(l),..., t{n — 1)) ; 

(4) there is an ‘operator’ F , acting on (possibly inhnite) sequences of power 
series, such that the recurrence relation for the t{n) ‘can be expressed 
in terms of generating series’, for example T(z) = r(T(z), T(z^),...) , 
which is abbreviated to T( 2 ;) = r{T(z)} . 

This is the most penetrating presentation we have seen of a foundation for recursively dehned 
classes of trees, a goal that we hnd most fascinating since to prove global results one needs 
a global setting. However their conditions have limitations that we want to point out. 

(1) requires t(l) = 1, so there is only one object of size 1; this means multicolored trees 
are ruled out. (2) indicates that one is using a specihcation^^ like T = {•} U •/Seq[y[(T) . 
The recurrence relation in (3) is the one item that seems to be appropriately general. It 
corresponds to what we call ‘retro’. (4) is too vague; after all, a function of (T( 2 ;), T( 2 ;^),...) 
is really just a function of T(; 2 ) since T( 2 ) completely determines all the T{z^). This 
formulation is surely motivated by the desire to include the MSet construction; perhaps 
the authors were thinking of ‘natural’ functions of these arguments like 
perhaps something of an effective nature, an algorithm. 

After this general discussion, without any attempt to prove theorems in this context, they 
turn the focus to simple classes T of rooted trees, namely simple classes are those for which 

(Ml) the generating series T( 2 ;) is dehned by a ‘simple’ recursion equation 

w = zA{w ), 

where A G M-°[[ 2 :]] with A(0) = 1. 

Additional conditions'^ are needed to prove their theorems, namely 

(M2) A{w) is analytic at 0, 

(M3) gcd{n e P : a{n) >0} = 1, 

(M4) A{w) is not a linear polynomial aw + b, and 
(M5) A{w) = wA'{w) has a positive solution y < Pa 

to guarantee that the methods of Polya apply to give the asymptotic form (★). (M2) makes 
T(z) analytic at 0, (M3) ensures that T( 2 :) is aperiodic, (M4) leads to < oo and T(/9 t) < 

have not needed a specification language so far in this paper—for this comment it is useful. Let • 
denote the tree with one node, and • /□ says to add a root to any forest in □. 

^^Their original 1978 conditions had a minor restriction, that a(l) > 0. That was soon replaced by the 
condition (M 3 ) —see for example |22| . 
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oo, and (M5) shows F(z,w) := w — zA{w) = 0 is holomorphic in a neighborhood of 
(pt, T(pt))- 

Thanks to the restriction to recursion identities based on simple operators they are able 
to employ the more powerful condition (M5) instead of our condition A(pa) = oo. Our 
condition is easier to use in practice, and it covers the two examples frequently cited by 
Meir and Moon, namely planar trees with A{w) = planar binary trees with 

A{w) = 1 + . For the simple recursion equations one can replace our A(pa) = cxd by 

the condition A'(pa) = oo. 


8. Open Problems 


(Ql) If T = E( 2 ;, T) with E G DOM[ 2 ;, tc] and T( 2 ;) has the shift periodic form 

then one can use the spectrum calculus to show there is an H G DOM[z, w] such that 
V( 2 ;) = H(z, V(z)). If E is open at (pt,T(pt)) does it follow that H is open at 
(pv, V(pv))? 

If so one would have an easy way of reducing the multi-singularity case of T to the 
unique singularity case of V. As mentioned in §12121 we were not able to prove this, 
but instead needed an additional hypothesis on E. Partly in order to avoid this extra 
hypothesis we used a detailed singularity analysis approach. 

(Q2) Determine whether or not the SetM operators can be adjoined to the standard operators 
used in this paper and still have the universal law hold. (See §ini21) 

A simple and interesting case to consider is that of identity (0,l,...,m)-trees with 
generating function T dehned hy w = z-f Does T satisfy the universal 

law? Example 78 shows the answer is yes for m = 2, but it seems the question is open 
for any m > 3. 

(Q3) Expand the theory to handle recursion equations w = G{z,w) with G having mixed 
sign coefficients. 

(Q4) Find large collections of classes satisfying the universal law (or any other law) that are 
recursively dehned by systems of equations 

wi = 

(**) : 

Wk = Qk{wi,... ,Wk), 
where the 0* are multivariate operators. 

Classes dehned by specihcations using the standard operators that correspond to 
such a system of equations are called constructible classes by Flajolet and Sedgewick; 
the asymptotics of the case that the operators are polynomial has been studied in 
ng, Chapter VII, provided the dependeney digraph has a single strong eomponent, and 
shown to satisfy the universal law (★). 

(Q5) The study of systems (**) is of particular interest to those investigating the behavior 
of monadic second order dehnable classes T since every such class is a hnite disjoint 
union of some of the 7) dehned by such a system (see Woods EDI). In brief. Woods 
proved: every MSO class of trees with the same radius as the whole class of trees 
satisfies the universal law. However his results seem to give very little for the MSO 
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classes of smaller radius beyond the fact that they have smaller radius. Here is a 
plausible direction: 

If T is a class of trees dehned by a MSO sentence, does it follow that T decomposes 
into hnitely many % such that each satishes a nice law on its support? 

(Q6) Among the MSO classes of trees perhaps the best known are the exclusion classes 
T = Excl(Ti,..., T„), dehned by saying that certain trees Ti,..., T„ are not to appear 
as induced subtrees. The Tj are called forbidden trees. A good example is ‘trees of 
height n\ dehned by excluding a chain of height n + 1; or unary-binary trees dehned 
by excluding the four-element tree of height 1. 

Even restricting one’s attention to the collection of classes T = Excl(T) dehned by 
excluding a single tree ohers considerable challenges to the development of a global 
theory of enumeration. 

Which of these classes are dehned by recursion? Which of these obey the universal 
law? Given two trees Ti,T 2 , which of Excl(Ti), Excl(T 2 ) has the greater radius (for its 
generating series)? 

From Schwenk IZ71 (1973) we know that if one excludes any limb from the class of 
trees then the radius of the resulting class is larger than what one started with. Much 
later, in 1997, Woods (201 rediscovered a part of Schwenk’s result in the context of 
logical limit laws; this can be used to quickly show that the class of free trees has a 
monadic second order 0-1 law. Aaron Tikuisis, an undergraduate at the University of 
Waterloo, has determined which Excl(T) have radius < 1. 

(Q7) Find a method to determine the asymptotics of a class U of free trees obtained from 
a recursively dehned class T of rooted trees. (See 11 17.H i 
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